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TOPOLOGICAL CURRENT 
IN FRACTIONAL CHERN INSULATORS 

TOHRU KOMA 


Abstract. We consider interacting fermions in a magnetic field on a 
two-dimensional lattice with the periodic boundary conditions. In or¬ 
der to measure the Hall current, we apply an electric potential with a 
compact support. Then, due to the Lorentz force, the Hall current ap¬ 
pears along the equipotential line. Introducing a local current operator 
at the edge of the potential, we derive the Hall conductance as a linear 
response coefficient. For a wide class of the models, we prove that if 
there exists a spectral gap above the degenerate ground state, then the 
Hall conductance of the ground state is fractionally quantized without 
averaging over the fluxes. This is an extension of the topological ar¬ 
gument for the integrally quantized Hall conductance in noninteracting 
fermion systems on lattices. 


1. Introduction 

The quantum Hall effect 031 El is one of the most astonishing phe¬ 
nomena in condensed matter physics. In fact, the quantization of the Hall 
conductance is surprisingly insensitive to disorder and interactions^ This 
robustness of the quantization reflects the topological nature of the Hall 
conductance formula. Namely, the topological structure never changes for 
continuously deforming disorder or interactions. 

The first discovery was that Thouless, Kohmoto, Nightingale and den Nijs 
[23] found that the Hall conductance is quantized to a nontrivial integer in a 
two-dimensional noninteracting electron system in a periodic potential and 
a magnetic held. Later, Kohmoto PI realized that the integral quantiza¬ 
tion of the Hall conductance is due to the topological nature of the Hall 
conductance formula. Namely, the integer is nothing but the Chern number 
which is the winding number for the quantum mechanical U( 1) phases of the 
wavefunctions on the magnetic Brillouin zone which can be identified with 
a two-dimensional torus. More precisely, the U( 1) phases are nontrivially 
twisted on the torus, and the winding number is given by the number of 
times when the phases rotate on the U{ 1) circle, and must be an integer. 

Although realistic systems have disorder and interactions, they treated 
translationally invariant noninteracting systems only. Instead of resorting 
to the topological argument in usual differential geometry, the method of 
noncommutative geometry [2J was found to be useful to show the integral 
quantization of the Hall conductance for noninteracting systems which do 
not necessarily require translation invariance [3j. (See also [21 EE IS IZO]-) 
However, its extension to interacting systems still remains an open problem. 


X A mathematical argument was given in m- 
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Another tricky approach which we will focus on the present paper was 
introduced by Niu, Thouless and Wu [25]. They considered a generic quan¬ 
tum Hall system which is allowed to have disorder and interactions. They 
imposed the twisted boundary conditions with angles 4 >i and 4*2 at the two 
boundaries in a two dimensional system, and assumed that the system ex¬ 
hibits a nonvanishing uniform spectral gap above a g-hold degenerate ground 
state. Under the assumptions, they showed that the Hall conductance aver¬ 
aged over the two angles is fractionally quantized. (See also [19].) 

Clearly, their method is artificial, and their result does not implies that 
the Hall conductance for fixed twisted angles is fractionally quantized. But 
one can expect that the twisted boundary conditions do not affect the quan¬ 
tization of the Hall conductance. Quite recently, Hastings and Michalakis 
[9] treated an interacting lattice electron system which does not necessar¬ 
ily require translation invariance, and showed that the Hall conductance is 
quantized to an integer irrespective of the twisted angles under the assump¬ 
tions that the system shows a nonvanishing spectral gap above the unique 
ground state. 

In the present paper, we consider interacting fermions in a magnetic field 
on a two-dimensional lattice with the periodic boundary conditions. As is 
well known, in generic lattice systems, a constant electric field is known 
to be useless in order to measure the conductance because the absolutely 
continuous spectrum of the unperturbed Hamiltonian often changes to a 
pure point spectrum. Therefore, one cannot expect that there appears an 
electric flow due to the constant electric field. 

In the present paper, we overcome the difficulty as follows: In order to 
measure the Hall current, we apply an electric potential with a compact 
support, instead of the linear electric potential or the time dependent vector 
potential which yield the constant electric field. Then, due to the Lorentz 
force, the Hall current appears along the equipotential line. Introducing 
a local current operator at the edge of the potential, we derive the Hall 
conductance as a linear response coefficient. For a wide class of the models, 
we prove that if there exists a spectral gap above the degenerate ground 
state, then the Hall conductance of the ground state is fractionally quantized 
without averaging over the fluxes. This is an extension of the topological 
argument for the integrally quantized Hall conductance in noninteracting 
fermion systems on lattices. 

The present paper is organized as follows. In Sec. El we describe the 
model, and state our main Theorem[l| In Sec. El we define current operators. 
Using these expressions, we derive the Hall conductance formula as a linear 
response coefficient in Sec. El Following the idea by Niu, Thouless and 
Wu, we twist the phases of the hopping amplitudes of the model, and we 
present Theorem El below about the stability of the spectral gap above the 
degenerate ground state against twisting the boundary conditions in Sec. El 
In Sec. the Hall conductance averaged over the phases is shown to be 
quantized to a fraction by using the standard topological argument. In 
Sec. El we consider deformation of the current operators, and prove that 
the statement of our main Theorem [T) is still valid for the deformed current 
operators so obtained. In other words, the fractional quantization of the Hall 
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condudctance is robust against such deformation of the current operators. 
Sections mm and m are devoted to the proof of the main theorem. The 
proof of Theorem [2] about the stability of the spectral gap is fairly lengthy, 
and therefore given in Appendix [Bj Appendices |A] and ICtill are devoted to 
technical estimates. 


2. Lattice Fermions in Two Dimensions 


Let us describe the model which we will treat in the present paper. Con¬ 
sider a rectangular box, 


r lWi 


r L ( 2 ) L (2)l 

2 ’ 2 

X 

2 ’ 2 


which is a finite subset of the two-dimensional square lattice Z 2 . Here both 
of L ^ and L ® are taken to be a positive even integer. We consider inter¬ 
acting fermions on the lattice A with the periodic boundary conditions. The 
Hamiltonian is given by 

( 2 . 1 ) Hq ^ = 'y ( tx,y c ]c c y + 'y ' y 1 U XliX2i __^ Xl n Xl n X2 ■ ■ ■ n Xl , 

X,y£A />1 Xi,X2,-..,Xi€:A 

where d,c x are, respectively, the creation and annihilation fermion opera¬ 
tors at the site x € A, the hopping amplitudes t XtV are complex numbers 
which satisfy the Hermitian conditions, 

t — /■* 

u y,x — u x,yi 

and the coupling constants U xliX2t ,„ jXI of the interactions are real; n x = dc x 
is the number operator of the fermion at the site x £ A. We assume that 
both of the hopping amplitudes and the the interactions are of finite range 
in the sense of the graph theoretic distance. 

Clearly, the Hamiltonian H^ of (12.11) commutes with the total number 
operator X^eA n z °f the fermions for a finite volume |A| < oo. We denote 

by Hq A ' ,N ' > the restriction of onto the eigenspace of the total number 
operator with the eigenvalue N. 

We require the existence of a “uniform gap” above the sector of the ground 
state of the Hamiltonian Hq A ' N ' > ■ Since we will take the infinite-volume limit 
A —» Z 2 , we keep the filling factor u = 1V/|A| to be a nonzero value in the 
limit. The precise definition of the “uniform gap” is: 


Definition 1. We say that there is a uniform gap above the sector of 
the ground state if the spectrum of the Hamiltonian Hq A,N ^ 

satisfies the following conditions: The ground state of the Hamiltonian 
Hq A N) is g-fold (quasi)degenerate in the sense that there are q eigenval¬ 
ues, E^\ ■ ■ ■ ,E { 0 N q \ in the sector of the ground state at the bottom of the 
spectrum of Hq A,N ^ such that 


E 


max 


){K 


(TV) 



as |A| —>• oo. 


Further the distance between the spectrum, 1 ,..., E^}, of the ground 
state and the rest of the spectrum is larger than a positive constant A E 
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which is independent of the volume |A|. Namely there is a spectral gap A E 
above the sector of the ground state. 


We prove: 


Theorem 1. We assume that the ground state of the Hamiltonian H^’ N ^ is 
q-fold (quasi)degenerate, and there is a uniform gap above the sector of the 
ground state. Further, we assume that, even when changing the boundary 
conditions, there is no other infinite-volume ground state except the infinite- 
volume ground states which are derived from the q ground-state vectors of 
the Hamiltonian Hq A ' N \ Then, the Hall conductance a 12 is fractionally 
quantized as 


( 2 . 2 ) 


1 p 

°V2 = -- 

2tt q 


with some integer p in the infinite-volume limit. 


The precise definition of the Hall conductance a 12 is given by (14.121) in 
Sec. [4] below. 


Remark . We can relax the assumption on the number of the ground states 
so that when twisting the phases of the hopping amplitudes at the bound¬ 
aries, there appears no other infinite-volume ground state at energies lower 
than or close to the energies of the q-fold ground state, except the infinite- 
volume ground states which are derived from the q ground-state vectors of 
the Hamiltonian Hq A ' N ' > . 

If one includes the charge e of the electron and the Planck constant h in 
the computations, then the Hall conductance (12.21) is written in the usual 
form, 


The possibility of the spectral gap above the q-fold degenerate ground 
state in quantum Hall systems was treated by m in a mathematical man¬ 
ner. The numerical evidences of the existence of nontrivial fractional Chern 
insulators were shown in [271 EH [261 ES]- 


3. Local Current Operators 

In order to measure the Hall current, we must define current operators by 
relying on the Ehrenfest’s theorem in quantum mechanics. More precisely, 
the velocity of the charged particle which gives the current can be determined 
by differentiating the expectation value of the center of mass with respect 
to time t. 

3.1. Current operators for a single particle. Consider the time evo¬ 
lution of a wavepacket of the single particle on the infinite-volume lattice 
Z 2 . The single-body Hamiltonian IKo corresponding to the present system 
is given by 


fK 0 if)(x) = t x>y ij;(y) + U x if(x) 




TOPOLOGICAL CURRENT IN FRACTIONAL CHERN INSULATORS 


5 


for if) € £ 2 (Z 2 ). The Schrodinger equation which determines the time evolu¬ 
tion of the wavepacket is given by 

at 

for the wavefunction at the time t. The expectation value of an observable 
a is 

_ ('•Pt, 
t ' (&,&) 

where the inner product is defined by 



x&? 


for tp, 'i/j € £ 2 (Z 2 ). 

The current due to the motion of the particle is determined by the velocity 
of the center of mass. The center of mass is nothing but the expectation 
value of the position operator. The position operator X = (X^ 1 ),^ 2 )) is 
defined by 

(X^ip)(x) = x^ip (x), j = 1,2, 

where we have written x = (x*' 1 ),^ 2 )) € Z 2 . The Ehrenfest’s theorem for 
the position operator is expressed as 

±(X^) t = (i[% 0 ,X ( %, 

where [A, B] is the commutator of two observables, A and B. Therefore, the 
natural definition of the current operator must be 

(3.1) j(i) = i[ft 0 ,X^}. 


Clearly, this is an infinite sum. We want to decompose the current operator 
into the sum of local current operators which are much more useful us to 
obtain the Hall conductance formula below. 

In order to define local current operators, we introduce an approximate 
position operator as 

(X^V)(:c) = f^ 3 \x)^{x), 
where the function on Z 2 is given by 


2 i 

fi\x) := k-l)-l 

k=i 


with the step function, 


9^\x;k) := 



x> k; 
x< k. 


Clearly, we have ||(X^) — X^)ip]\ -A 0 as t —>• oo for (X^'ip) € £ 2 (Z 2 ). 
By replacing Xk>) with X^ in the right-hand side of (13.11) . we obtain the 
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approximate current operator as 

3® =i[M 0 ,X? ) ] 

21 

= X>[K| o,0 U) O- ;*-*)]■ 

k =1 

The summand in the right-hand side is interpreted as the local current op¬ 
erator. Namely, the local current operator J^\k) across the k- th site in the 
j direction is given by the commutator of the Hamiltonian ‘Kq and the step 
function 9^\- ■ ■ ; k) as 

3.2. Current operators for many fermions. The step function is writ¬ 
ten 

(3.2) flW(fc) = 9^\x-k) n x 

x£l? 

in terms of the number operator n x of the fermions. Therefore, the local 
current operator J^\k) is formally written 

(3.3) J ij \k) = i[H^ 2 \9 {j \k)} 

= i y ' y ' tx,y[ c x c y, n z\ 
x,yeZ 2 z&? 

z^>k 

= i 'y ] y ] (txyCxCy — ty X CyC X ), 

x&? y&L 2 
x^<k y U)>k 

where we have used the expression m of the Hamiltonian, and we have 
written x = Although this is clearly an infinite sum, the ex¬ 

pression can be justified for a localized wavefunction of many fermions. For 
the present finite-volume lattice A, we can define the local current operator 
jW(k) as 

{k) := i y ^ y ^ ( t X ycl.Cy ty X c!yC x ) 

yGA 

x^<k yU)>k 

with the periodic boundary conditions because the range of the hopping 
amplitudes is finite. In order to measure the Hall current, we introduce an 
approximate local current operator as 

(3.4) 4 1} (M) 

: = i 'y^ 'y^ 'y^ 'y^ ( txyc x c y — t yx c^c x ), 

xW<k t-Ji<xV 2 )<e.+'H yW>k e-N<y( 2 ) <e+N 

where N is a positive integer and £ is an integer. Namely, N is the cutoff in 
the second direction. 
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4. Linear Response 

In order to measure the Hall current, we apply an electric potential to the 
present system. For this purpose, we introduce a 2M x 2M square box as 

r u(k,£) ■= {z = (z (1) ,z (2) ) \k-M< X W < k + M,£< x ( ^ <£ + 2M} 

for a positive integer M, and the sum of the number operators n x on the 
box as 

(4.1) x(r m(M)) : = Hx - 

xer M (k,e) 

The latter yields the unit voltage difference from the outside of the region. 
The current operator segment J^\k,£) of (13.4h goes across the bottom side 
of the square box Tj^i(k, £). When we change the potential energy inside the 
box by using the potential operator x(Tjvt(fc, £)), the Hall current is expected 
to appear along the boundary of the box Tj^(k,£) due to the Lorentz force. 
The Hall current is measured by the observable J^\k,£). More precisely, 
we use a time-dependent electric potential. The Hamiltonian having the 
electric potential on the region T j^(k,£) is given by 

(4.2) H^ A \t) := H { 0 A) + XW{t) 
with the perturbed Hamiltonian, 

(4.3) w(f) = e ^ x (r M (M)). 

Here the voltage difference A is a real parameter, and the adiabatic param¬ 
eter ij is a small positive number. We switch on the electric potential at the 
initial time t = —T with a large positive T, and measure the Hall current 
at the final time t = 0. 

The time-dependent Schrodinger equation is given by 

= H^ A \t)^ N \t) 
dt 

for the wavefunction 4 >( N \t) for the N fermions. We denote the time evo¬ 
lution operator for the unperturbed Hamiltonian Hq A ^ by 

(4.4) U^ A \t, s) := ex.p[—i(t — s)Hq A ^] for t,sG R. 

We choose the initial vector $( 7V )(—T) at t = — T as 

^(iv)(_ T ) = o)$W 

with a vector Then, the final vector ^ N ) = = 0) is obtained 

as 

/ o 

ds C/q A) (0, s)W (s)Rq A) (s, 0)$^ + o(A) 

by using a perturbation theory m, where o(A) denotes a vector 4/^ with 

the norm ||4'^|| satisfying ||4/^||/A -A 0 as A -A 0. 

We denote the q ground-state vectors of the unperturbed Hamiltonian 
H {A ’ N) by with the energy eigenvalue m = 1, 2,..., q. We choose 

the initial vector 4?^ = 4 >q^ with the norm 1, and write the corresponding 
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final vector at t = 0 as 'I'W = Then, the ground-state expectation 

value of the approximate local current operator J^\k,£) is given by 


(4.6) 




m= 1 


Let be the projection onto the sector of the ground state of Hq A,N ^ , 


and we denote the ground-state expectation by 


w. 


(A ,N) 


(-):=;Tr(-)n 


(A ,N) 


where q is the degeneracy of the ground state. Using the linear perturbation 
(14.51) . the expectation value of (14.61) is decomposed into three parts as 


where 


<4 1} (M)> = (J$\k,e)) 0 + {j£\k,l)) 1 + o(A), 

(j£\k,e)) 0 = 4 A ’ N \j£\k,e)), 


and 

(4.7) (J$\k,£)) 1 

rO 

= -i\ I dsu {A ’ N \j^\k,£)U^ A \o,s)W{s)U^ A \s,0)) + c.c. 


The first term (j^\k,£)) 0 is the persistent current which is usually van¬ 
ishing. We are interested in the second term ^ ,7^ (k, ^)) 1 which gives the 
linear response coefficient, i.e., the Hall conductance. We write 

(4.8) X (A) (T m (M); S ) := ^ A) (0, s) X (r M (k,£))U^ (s, 0). 

Using this and the definition (14.31) of W(t), the contribution (14.71) of the 
expectation value of the local current is written 

(4.9) (j£\k,£)) 1 =i\ J° T dse* s u;^ N \[x {A \r M (kJ);s),j£\k,£)}). 

Note that, by using integral by parts, we have 
r° 

(4.10) 


' —T 


ds e vs u { 0 A ' N \[x iA} {T M {k,£);s),J^\k,£)}) 
se^J A ’ N) ([ X (A) (T M (k, £)■ s), 4 1} (k, £)]) 
ds ?7se r?s 4 A ’ Ar) ([x (A) (r^vc (k, £);s), (k, £)]) 


- o 

. -T 


i: 


' —T 


ds se vs uj o A,JV ' ) 


^-x (A) (r M (k,£y,s),jW(k,£) 


ds 


Clearly, the first term in the right-hand side is vanishing as T -A oo. The 
second term is also vanishing as r] —» 0 after taking the limit T -A oo. We 
show this in Appendix 0 The third term in the right-hand side is written 

^ ds seV s J A ’ N \[j£\k,£),J^(T M (k,£y,s)}), 


l-T 
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where 

j^ A \r M (k,£)-s) 

■= ^x (A) (r M (M);s) = u^\o,s)i[H { 0 A) ,x(,r M (k,mu^\s,o). 

For getting the second equality in the right-hand side, we have used the 
definition (14.4p of the time evolution operator U^\t,s) and the defini¬ 
tion (14.81) of the operator x^Ogm (k,£)‘,s). Clearly, the above operator 
J l ' A \Tj^(k,£)-, s) is interpreted as the current across the boundary of the 
region Tj^(k,£) at the time s. From these observations, we define the Hall 
conductance for the finite lattice A as 

(4.11) a[ A ’ N \r,,T,X, M) 

:=i /° ds se 713 ^’^ ([J^ (k, £), J (A) (T M (k, £); s)]) 
J-T 

because A is the voltage difference. The Hall conductance in the infinite- 
volume limit is given by 

(4.12) a 12 ■= lim lim lim lim lim T, N, M). 

N —>oo M—>oo ? 7 —s-0 T—voo A /'l? 


5. Twisting the Phases 


Consider 

Hf \(f>j,k) := exp[— i(f>jd^\k)]H^‘ ^ exp [i<j)j6^\k)\, for cf>j £i, 

where 6^\k) is the step function of m- This transformation changes the 
hopping amplitudes as 

t xy —*■ t xy e^ j for < k, > k, 

and the opposite hopping amplitude t yx is determined by the Hermitian 
condition tyx ~ t% y - The rest of the hopping amplitudes do not change. In 
both of the first and the second directions, we define the twist by 

H { 0 Z 5 (0i, fci; 02, fc 2 ) 

:= exp[—i((/>i0^ 1) (fci) + H { 0 Z ) exp^^id^^ki) + ■ 

This transformation is justified for localized wavefunctions. 

Relying on this rule, we can twist the phases of the hopping amplitudes 
of the finite-volume Hamiltonian H^ A \ We denote by H^- 1 (r/>i, k \; fa, ^ 2 ) 
the corresponding finite-volume Hamiltonian on A. When both of k\ and 
^2 are placed at the boundaries of the lattice A, the boundary conditions 
of the system become the usual twisted boundary conditions. We write 
Hq A \(/) i;<fe) f° r the Hamiltonian with the twisted boundary conditions for 
short. Further, one has 


H (f - 1) = exp 

^ ^ Tl x 

Hq 1 \(/>j,k)exp 

■ " 

nx 


xez 2 



xez 2 


x^=k —1 2 (-4 =k—\ 
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Thus, one can change the position of the twisted hopping amplitudes by 
using the unitary transformation which is local in the 7 -t.l 1 direction. 
Similarly, we can consider 

fci; 0 2 , k 2 ) 

:= exp [-i (£q) + (fc 2 ))] (k) exp [*(fci) + (fc 2 ))] 

for the local current operator J^\k). Therefore, in the same way as in the 
case of the Hamiltonian, one can change the position of the twisted phases 
as 

\,ki - 1 ; 02 , ^ 2 ) 


= exp 

8 

g 

W 

* <S> 

1 

fci; 0 2 ,fc 2 )exp 

w 

1 * ^ 


ccEZ 2 




x^=ki— 1 


=k\ — 1 

Further, one 

has 




jnr H o‘ 2 \ < l > j 'k) = ^p[-ifijO {3 \k)]i[H^ 2 \e^\k)]eyi.p[ifi j e^\k)\ 

Ofpj 

= exp[— i(j)j6^\k)\ J^\k) exp [i<f>jd^\k)\. 

Namely, we can obtain the local current by differentiating the Hamiltonian. 
Thus, one has 

r\ 

(5.1) —H^fa, kr,(hM) = J < ' J \k j -,(t>i,k 1 -,(j)2,k2), j = 1,2, 
for the finite-volume Hamiltonian Hq A \ 

Theorem 2. Under the same assumptions as in Theorem\T) we have that 
the ground state of the Hamiltonian Hq A \(/)i, fa, & 2 ) having the twisted 
phases <f\ and 2 at the positions k\ and L 2 in the first and the second di¬ 
rections, respectively, has the same q-fold (quasi)degeneracy as that of the 
Hamiltonian H^ without the twisted phases. Further, there exists a uni¬ 
form spectral gap above the sector of the ground state for any <fi and <f 2 , and 
there exists a positive lower bound of the spectral gap such that the bound is 
independent of the phases 4 >i, 4 > 2 - 

The proof is given in Appendix [Bj 

Remark . As mentioned in Sec. [TJ the method of twisting the phases at 
the boundaries was first introduced by [25] for the quantum Hall systems. 
Later, the method was extended to strongly correlated quantum systems 
such as quantum spin systems in mmm, in order to measure the local 
topological properties such as a local singlet pair of two spins. 

6 . A Topological Invariant 

Consider the Hamiltonian Hq A \()i, k] </> 2 ,^). Relying on Theorem O we 
denote the q vectors of the ground state by $ 2 , 2 ) with the energy 

eigenvalue E^f^ ( 0 i; </> 2 ), m = 1 , 2 ,..., q, and denote the excited-state vectors 
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, k; cf> 2 ,i) with the energy eigenvalue 4 ^( 01 ; fa), n > 1 . We write 
0 = (0i,02)? and 0 = ((j>i,k\(j) 2 ,£) for short. We define 


( 6 . 1 ) a[ 2 ,N \&) 


EE 


^ (1) (fe; 0 )^^ ( 0 )) ( 0 ), ^ (2) (^; 


m=ln>l L (4^(0)-^ iV) (0)) 2 

_ (4^(0), J( 2 ) (l; 0)44 (0)) (44 (0), (fc; 0)<Q0)) 

(4^(0)- 4*° ( 0)) 2 

We write 42(0) for the vectors of the ground state of the Hamiltonian 
44 (0U 02 ) with the twisted boundary conditions with the angles 0j and 
02, and 44(0) for the vectors of the excited states. We define 


( 6 . 2 ) 42 -A,, (<W 


-EE, 


(42(0), JW (fc)44 (0)) <44 (0), J (2) W4^(0)) 


0,m' 


(42(0)-^ JV} (0)) 2 

_ <4^(0)^ (2 H^)4 jV) (0))(4 JV) (0),j^(fc)4!S(0)) ' 

( 42 ( 0 )- 4 *° ( 0)) 2 

This is the standard form of the Hall conductance for the degenerate ground 
state with the twisted boundary conditions. Since we can change the posi¬ 
tion of the twisted hopping amplitudes by using the unitary transformation 
as shown in the preceding section, this conductance 012 ’^ (0) * s e( l ua l to 
the above 4-2’^(0) °f (EH)- 

Niu, Thouless and Wu )25j showed the following: When averaging the 
conductance cr^’ N \^>) over the phases 0 1 and 02, the averaged Hall con¬ 
ductance is fractionally quantized as in (16.121) below. Along the lines of [T9], 
we shall give the proof. 


Using a contour integral, the projection Pq 
ground state is written 


(A,TV) 


(0) onto the sector of the 


P 


(A ,1V) 


( 0 ) = 


1 

2iri 


dz- 


z-H, 


(A,TV) 


( 0 ) 


because the spectral gap exists above the sector of the ground state as we 
proved in Theorem [2] Note that 


P, 


(A,TV) 


0,1 


( 0 ) : = 


d 


P 


30T° 

1 

27 n 

1 


(A,TV) 


( 0 ) 


dz- 


z — H, 


(A,TV) 


(0) L 


dH { 0 A ’ N \& 

50i 


z-H, 


(A,TV) 


( 0 ) 


= x—: 0 dz 


2mJ z -H (A ' N) (< t>) 


«/ ( 1 ) (^; 0 )- 


z- H, 


(A,TV) 


( 0 )' 
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where we have used (15.ID . Similarly, one has 



Using these, we have 



Remark . In the infinite-volume limit, A Z 2 and N /* oo, one formally 
has 



where we have used 



which are derived from (13.31) . This expression of the Hall conductance 
djf ’°°' > has the same form as that for formally applying the method of 
noncommutative geometry mmm to an interacting fermion system. Ac¬ 
tually, the expression can be obtained by replacing the single-fermion step 
functioifl with that for many fermions. But, justifying the expression in a 
mathematical rigorous manner remains a challenge because the many-body 
step function 6^\kj) is expressed in terms of the infinite sum of the number 
operators. 

The following proposition is essentially due to Kato m- See also Propo¬ 
sition 5.1 in [El- 

Proposition 3. There exist orthonormal vectors m = 1 , 2 ,...,^ 

such that the sector of the ground state of ((f) is spanned by the q 

vectors (&)> ani ^ th a t a H the vectors are infinitely differentiable 

with respect to the phase parameters <f> £ [0 , 27t] x [0, 2n\ . 

Relying on this proposition, the right-hand side of (16.31) is written 



m =1 



>(A,iV) 

0,2 




m= 1 



2 See, e.g., [21I51I2U]. 
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except for the factor i/q. Note that 


>(A,IV)/ 


= jst'wCw+ 


(A,IV), A 3 i (IV), 


d< t>3 


for j = 1, 2. Therefore, one has 


(6.4) 


1 -P O (A>JV) ( 0 ) 


0 


C!(0)= f ° r j = 2 


,(A.JV) 


(IV), 


J ab¬ 

using these identities, the above trace is written 


TV P r 


(A,IV) 


( 0 ) 


d(A,JV) 

0,1 


m=l 


(0),^ ,JV) (0) 

d 


Ed- ^(*)]^*g2(*), [i - 


00 i 


002 


- <[1 - ^(*)l £*&!(*). [1 - pf 

-£{(£*»! 

m =1 ^ ' 

5 n ofA.Jvj/^M 5 c(jv) 


0 


00 2 


iEw.[i-Po (A ' N) w]g^*Ew 


E 


ra=l 


9 w)fi\ 9 *m 




, 00 ! U ’ mV " /5 002 
where we have used the identity, 

i(A) 




= dd+(dw. ^-dw 


for j = 1 , 2 , and m,m' = 1 , 2 , ...,</, in order to get the third equality. 
Consequently, the Hall conductance ’ ,N ^ (0) of (16.31) is written [19] 


(6.5) a 


^(0) 


0 


0 


0 


0 


- ; E ^)> - 

Next, following m, we show that the Hall conductance averaged over the 
phase 0 is fractionally quantized. We define the averaged Hall conductance 
as 

1 


( 6 . 6 ) 


a 


(A,IV) 
12 


( 0 ) := 


(2vr) 


0 M 0 2 


' [0,2tt] X [0,2tt] 

Using the expression (16.51) of the Hall conductance, the integral in the right- 
hand side is written 


(6.7) 


[0,2ir] X [0,27 t] 


d0ld02 <7 


( A’IV)( 0 ) = * /(l)_/( 2 ) 


12 


9 L 
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where 


r2n 

/(1) = E / 


m=l 


I 2 *, <fc) ■ I 2 *, <fc) 




and 

/ (2) = £ / #1 


^ /*27T 


m=l ' 7 ° 


®SS (* . 2,r ). .a 


-(i<5(*,o),A$W(*,o) 

Here, we have dropped the fc, ^-dependence of the vectors $om(^i> 02, 

for short. Clearly, the set of the ground-state vectors, ^omWh.’27 t), m = 
1,2 ,... , g, connects with that of 0), m = 1,2 ,..., q, through a q x q 

unitary matrix C^ 2 \(j)i) as 


(6-8) ^oX(0 1 > 27r ) = E for m = 


d 2 ) 


(TV) 


m'=l 


where ) are the matrix elements of the unitary matrix C^ 2 \(/)i) 

By differentiating this, one has 

= a r»ci 2 L.w,) 

m'=1 . 


3</>i 


*©(*>.») + 


for m = 1,2,..., q. By using these, one obtains 

E (^o5(^l ; 27r )>^ $ 0,m(^l’ 27r )i 


ra=l 


q q q _ 

EE E 


C (2) ,(<fc 

m,m' 


„aci 2 ,>„„(*) /A(K) 


m=l ?n'=l m"=l 
-<( 2 ) l A. 


^^<CJAi,o), <„!»(*, o)> 


E E c ■ 


sc< 2 > ,(*,) 


'=1 m=l 


3</>! \ 0,m 

Substituting this into the expression of , we have 

/ {2) = / 2 "#i Tr 

Jo d( Pi 
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Since the unitary matrix i) can be expressed in terms of a Hermitian 

matrix 0^ 2 ^((/>i) as 

C (2) Oi) = exp [z0 (2) Oi)], 


we 


obtair@ 


f2w c\ 

(6.9) lW=i #1 Tr — 0( 2) (<^i) = irTr0( 2) (27r)-Tr0( 2 )(O)l. 

Jo OtJ) 1 

Similarly, one has 


<7 



m '=1 


with a unitary matrix C^\(/> 2 ). Therefore, in the same way, we obtain 
/ (1 ) = i [Tr 0( 1) (2vr) - Tr 0^(0)], 
where we have written 


CWfa) = exp [z©^ 1 ^ (<^>2)] 

in terms of the Hermitian matrix ©^'((^l). From (16.611 and (16.711 . we obtain 
[Tr 0( 2 )(2tt) - Tr 0^(0) - Tr 0 ( 1 )(2vr) + Tr © ( 1 ) (0)1. 

(27 T)-q 

Therefore, it is sufficient to show 

(6.11) Tr 0 ( 2 ) (2t r) - Tr 0 ( 2 ) (O) - Tr 0 ( 1 ) (2tt) + Tr 0 ( 1 ) (O) = -2vrp 

with some integer p. 

For this purpose, we denote by 6 q V ^(^i, <j > 2 ) the g-component vector with 

the m-th component, 6 [/^(</>i, 62 ), m = 1,2,..., q. Then, the relation (16.811 
is written 

^ 2tt) =C( 2 \cl> 1 )& 0 N) (</>!, 0 ). 

Similarly, the relation (I6.10H is written 

df ) (27r> 2 ) = cW(0 2 )4>^ ) (o,< ? i 2 ). 

For (f>i = 02 = 27r, one has 

6^(271,271) = C (2) (271)6^(271,0) = C (1) (271)6^(0,271). 

For 4>i = 4> 2 = 0, 

6^(0,271) =cW{Q)*f> (0,0) 

and 

6^(271,0) = c ( W(0)6f ) (0,0). 

From these equations, the following relation must hold: 

C (2) (2ti)C (1) (0) = C (1) (2ti)C ,(2) (0). 


Taking the determinant for both sides of this equation and using the formula 
detexp[i0( J )(0)] = exp[iTr 0^ (</>)] for j = 1,2, one has (16.111) . 
Consequently, one has 


( 6 . 12 ) 


a 


(A,TV) 
12 


( 6 ) := 


(2ny 


K 


[0,2vr] X [0,27r] 




( 6 ) 


_L p 

271 q 


^The differentiability of 0^(0i) is justified in Appendix [Hj 
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with an integer p. Since a[ A ’ N \(f>) = d | AN \&) as mentioned above, we 
obtain 


(6.13) 


1 


[0,27r]x[0,27r] 


# 1#2 ^2 


1 p 

2i r q 


for the Hall conductance a^’ N \<f)) with the simple twisted boundary con¬ 
ditions. 


7. Topological Currents 

We write H^ A \cj)) for the Hamiltonian Hq A \c/)i, k\ 4>2,£) in Sec. [5j Con¬ 
sider a further deformation of the Hamiltonian H^ A \^>) as 

HQ A \cj),a) := exp[-ia0in a; ]i7g A) (^)exp[ia^in x ], 

where we have written a = (a, x) for short. Then, the relation (ED is 
modified as 




= exp[-ia(j)in x } { J (1) (fc; 0) + ia[H { 0 A \&), ?i x \j exp^a^n^ 


We write J 9 \k; <f>,a) for this right-hand side. Namely, we have 
(7.1) -^- H { A \ j >, a ) = jU ( hA , a ). 

Clearly, this operator (k; <fi, a) is the local current operator which is 
obtained by twisting the phases of 


(7.2) (k, a) := J (1) (k) + ia[H {A) , n x ] 


with the angles (j>i and <p<i- In the infinite volume, the operator is formally 
written 


J^\k,a) = i[Hf\e^\k)} + ia[Hf\n x } 
= i[H^ \o 9 \k) + an x \. 


That is to say, the step potential is slightly deformed with a at the site x. 

In the same way, the projection operator onto the sector of the ground 
state is transformed as 

PQ A ' N \(/),a) := exp[-fa^in a: ]Pd A ’ Ar) (0)exp[iQ!</>in a: ] 

1 / dz 

2m J z - H {A ’ N) (j,,a) 


We write 


and 


Then we have: 


P$ N \<l>,a) := 


9 d(A,JV) 


d<h 0 


^ 0,2 W,a)- d(h r Q 


(&,&) 

(0,d). 






TOPOLOGICAL CURRENT IN FRACTIONAL CHERN INSULATORS 


17 


Lemma 4. 


[0,2ti-] x[0,2tt] 


d(j)id(j )2 Tr PQ A ’ N \(j),a ) 


Pc 


(A,IV) 


[0,27t] x [0,27r] 


# id<h Tr P O (A>JV) (0) 


Pc 


0,1 

(A, IV) 


0,1 


{^,a),P^ N] (^,a) 


The proof is given in AppendixlH This lemma implies that the topological 
invariant, the Chern number, does not change for the local deformation 
of the phases of the Hamiltonian. But the local current operator in the 
first direction changes from (13.3D to (17.21) . Even when starting from this 
deformed local current operator J^){k,a), our argument holds in the same 
way. In consequence, we can prove the same statement as in Theorem [T] for 
the deformed local current operator. Since we can repeatedly apply local 
deformations, we can define more generic local currents in the following way. 
Consider generic functions d^\x m ,k,£), j = 1,2, which satisfy 

£} < ' 1 \x-,k,£) = 6^ 1 \x-,k), d^ 2 \x-,k,£) = 6 <y2 \x,£) for dist(x, (k, £)) > Ro 

with a large positive constant Rq. Namely, the functions , d^\x' 1 k,£) coin¬ 
cide with the step functions at the large distances. At short distances, the 
functions •d^\x]k,£) can take any real values. We define 


•d^\k,£) := ^2 d^\x',k,£) n x . 

x&? 


Then, the twisted Hamiltonian is formally given by 

\<f>) :=exp[- } exp[^ (k,£)\. 

3 = 1,2 7 = 1,2 

The generic local current operators are 

3 U) (k,£) :=i[Hf\d^\k,£)]. 

The potential operator of dUB) can be also replaced by a generic electric 
potential such that, except the boundary ribbon region, the potential takes 
the unit inside the finite large region and zero outside region. Namely, the 
potential coincides with the characteristic function of the finite region at the 
large distances from the boundary of the region, while the potential takes 
any values at the short distances from the boundary of the region. 


8. Estimates of the Current-Current Correlations 

The conductance (14. lip is written 

/ o 

ds se vs 


x - £ £ [<C, 


m= 1 n> 1 


- ($w,j(t m (a : ^))^ jv ))(^ jv ),4 1) ( 

in terms of the ground-state vectors with the eigenvalue Eq^, m = 
1,2,... ,q, and the excited-state vectors with the eigenvalue E^\ n > 









18 


T. KOMA 


1. Here, we have dropped the twisted phase dependence for short. Consider 
the limit 

^ 12 ^ (0) °o, N, JVC) := lim lim T, N, M). 

?7—s-0 T—too 

Using the integral formula (lA.lj) in Appendix [A] this is computed as 
4^(0, oo, Af, M) 


EE 

m= 1 n> 1 


(4 N 2-En N) y 

,{N) _ j?(N)\2 


\^0 ,m 

where we have written 

j(r M (M)) := J (A) (r M (M); 0 ) 

for short. In the same way as in Appendix [Aj we have 


(8.1) cr^^T.N.M) -cr&’^oo^M) 


||4 1) (fe,£)|| ||j(r M (fe,£))||. 

Both of the current operators, J^\k,£) and J(Tj^(k,£)), are written as 
a sum of local operators. We want to show that the dominant contributions 
to the conductance in the double sum are given by local operators in the 
neighborhood of the point (k,£). In other words, if the distance between 
two local operators in the sums is sufficiently large, then the corresponding 
contribution are negligible. 

We define 


< 2 


2A E + r? 1 + TAE 


AE 4 


-rj + 


A E 2 


~VT 


1 


«AB»:= AEC."? 1 ) 


q 

9=1 n> 1 

for local observables A and B. 


( _ p (-^0 \2 

VUO.m ) 




Proposition 5. The following bound holds: 

I i(A', B))\ < Ce~ Kr , 

where r = dist(supp A, supp B), and C and k are positive constants. 
Proof. We define 

B := B - P^BpW) 

and 

B^\z) := e iH o A)z Be~ iH o A)z for z € C. 

Using the integral identity 


1 

A £2 


f 

J o 


ds s e 


—A£s 
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with A£ > 0, we have 


«A B » 

1 ^ /*oo 

:= \ E E / 

i ^ rcr 

= \ E E / s(^.M JV) )exp[-(4 JV) 

®m=ln>l 



ds s(<h 


(IV) 
0 ,m> 


yKff’jexpl-tEW 




where c is a positive constant. Clearly, the second sum leads to the desired 
bound because of the spectral gap above the ground state. The first sum is 
written 

ds s Uq A,N \AB ( ' A ' > (is)). 

Therefore, it is sufficient to estimate uiQ i ’ N \AB(ib )) for 0 < b < cr. By 
definition, one has 

J Q A ’ N \AB^ A \ib)) = J A ’ N) (AB^ A \ib)) - u iA ’ N \AP^ N) B^ A \ib)) 

where we have written 

B^\z) := e iH o A)z Be~ iH ° A)z for z £ C. 

For a large distance r = dist(supp A , supp B), we can expect the exponential 
clustering 

J A ’ N) (AB^ A \ib)) ~ J A ’ N) {AP^ N) B^ A \ib)) + 0(e~ Kr ) 

for a small b because of the spectral gap above the ground state. Actually, we 
can prove this by using the Lieb-Robinson bounds in the same way as that 
for evaluating the integrand in (1B.24[) in Appendix lB.21 (See also Theorem 2 
in [23].) □ 



Using Proposition [5] we have 


012 "' N) (0, °o, N, 0) ~ <^ 2 (0) 


= 0 


(8.2) lim lim lim 

N—>-oo M —>-00 A /'I? 

for the Hall conductance cr[ A,N ' > (0) of (16.21) with the twisted boundary phase 
0 - 


9. Twisted Phase Dependence of the Hall Conductance 
We write 

a[ A ’ N \r,,T,X,M,<j>) 

for the conductance of (14.111) in the case of the twisted boundary conditions. 
We write 


A^ A \t) \= e ltH ° A) Ae ltH ° A> 










20 


T. KOMA 


for the time evolution for a local operator A. Let 12 be a subset of A, and 
write Hq^ for the Hamiltonian restricted to the subset 17. Let A be a 
local observable with supp A C 12, and define 

A™(t) := e itH ° n) Ae- itH ^. 

This is the time evolution of A on the region 17. Then, one has 

(9.1) ||A A, (t)- A n >(t)|| < /‘*||[(4 A, -A n, ).A n >(t-»)]|| 

JO 

for t > 0. When the distance between the supports of two observables 
— Hq^) and A is sufficiently large, the norm of the commutator in 
the integrand in the right-hand side becomes very small for a finite t. This 
inequality can be proved by using the Lieb-Robinson bounds in Appendix ICl 
We define 

(9.2) 0 ^ 2 ,iV) ( 77 , T, N, M, 17, 0) 

:= i f ds S e^ s 4 A ’ N) ([jU ( k , e), (r M (fc, ey «)] ■ 0 ). 

This is given by replacing the current operator s) by 

J (q) (T m (M);s) := e itH ^ J(r M (k,£))e~ itH o n) 

in the expression (14.1111 of the conductance ' N ^ ( 77 , T, AT, M, 0) in the case 
of the twisted boundary conditions. Note that 

(9.3) d[ A ’ N) ( 77 , T, N, M, 0) - a[ A ’ N) ( 77 , T, N, M, 17,0) 

= * dsse** 4 A,JV) ([4 1} (*> 4 J(A) eys)- j w (r M (k, ey «)] ; 0 ). 

J-T 

Using the above bound (EH), we can prove that the difference between the 
two conductances becomes very small for a large lattice A. Actually, we 
choose the region 17 so that 

dist(17, dA) = 0 (L) 

and 

dist(ai7,T M (fc,£)) = 0 (A), 

where dA and 317 denote the boundary of the regions A and 17, respectively, 
and L = minjLW, L^}. Then, we have 

(9.4) of^’ JV) (77,T,N,M,0) -d^ ,iV) (77,T,N,M,17,0) -> 0 as A /> Z 2 

by using the Lieb-Robinson bounds (IB.121) below in Appendix lB.il (See also 
[2H 013(23].) 

In the right-hand side of (19.21) . the support of the commutator of the two 
current operators is a finite subset of A and apart from the boundary of A. 
Therefore, we can expect that the effect of the twisted boundary condition 
is exponentially small in the size L in the approximated Hall conductance 
(19.21) . Actually, we can prove the inequality (19.51) below. 
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To begin with, we note that the expectation value of a local operator A for 
the ground state of the Hamiltonian H^ A ' N \cf)) with the twisted boundary 
conditions is written 


4 A ’ iV) (H;0) = iTrHP o (A ’ Ar) (^), 

q 


where PQ A,N \cf)) is the projection onto the sector of the ground state. Using 
contour integral, one has 

1 


P, 


(A, IV) 


<« = 5S / 


dz- 


z-H, 




Further, by relying on the existence of the spectral gap above the sector of 
the ground-state, one obtains 

1 „ , 1 


9 p(A,JV) 

JP n 


9 h 

where we have written 


<« = * hj 


dz- 




Bj (<!>)- 


z-H, 


(A, IV) 


( 0 ) 




From this, we have 


-,^-Tr d? 0 (A ’ A) ^\ - — 

9( t>j 


dz TV A- 


z-H, 


(A, IV) 


( 0 ) 


-BMY 


z-H, 


(A,IV) 


(0) 


Integrating both side in the case of j = 1, we obtain 


Tr AP^icj) i, &) - TV AP 0 ^’ JVJ (0, 0 2 ) 

1 /' 1 


,(A,IV). 


f 


dz Tr A- 




z-H, 


(A, IV) 




2 ™/ " ^z - H^' N \(t>' 

where we have written <fi' = Therefore, 

u ( 0 A ’ N) {A; &) - w^ A,Ar) (H;0,(/>2) 


/•(pi -i ^ ___ 

[ #'-V E 




452 wo 


+ (H -B- Bi(<t>')). 


This right-hand side can be evaluated in the same way as in the preceding 
Sec.0 Thus, the difference of the two conductances with the different phases 
is exponentially small in the linear size L = minjL^, L^} of the lattice A 
as 


^ (r?, T, IT, M, n, d>i, fa) - cr£’ JV ' ( 77 , T, N, M, SI, 0,0 2 ) 


(A, IV), 


where 


e(rj,T) : = 


< C(r/, T) x 0(exp[—Const.L]), 


ds se vs 


' —T 
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Iii the same way, we obtain 


(9.5) 


^2 ’ N) (v, T, n, m, n, fa , fa) - <x£’ JV J (n, t, n, m, n, o, o) 


(A,iV), 


< ^2 N) (v,t, n,m, a, fa,fa) - a£’ N) (r],T,X,M,Sl,0,fa) 

+ \a £’ N) (77 , T, N, M, 0,0, 0 2 ) - ef£ )JV) ( 77 , T, X, M, ft, 0,0) 

< C( 77 ,T) x 0 (exp[—Const. L]). 

This is the desired inequality. 


10. Proof of Theorem Q] 

Now, we shall complete the proof of Theorem [1] The difference between 
the two conductances with the twisted phases, 0 = {fa, fa) and ( 0 , 0 ), is 
evaluated as 


(10.1) 

N) (v, T, N, M, fa, fa) — d£ N) ( 77 , T, N, M, 0,0) 


< 

N) 0 V , T, N, M, fa , fa) - d£ N) ( 77 , T, l N,M,Q,fa,fa) 


+ 

(v,T,X,M,n,fa,fa)- ~^ N) ( 77 ,T, N, M, ft,0,0) 


+ 

d[2 ,N) ( 77 , T, N, M, ft, 0,0) - d[2 ,N) ( 77 , T, N, M, 0,0) . 


From (19.41) and (19.51) . the right-hand side is vanishing in the limit A /* 7?. 
By using the formula (16.131) for the averaged Hall conductance, we have 


a 


iJ’^r^M^o) —^ 

Ztt q 


( 27r ) 2 l[ 0 , 


2i r) X [0,27r) 


d(j)\d<j >2 ^[ 2 ’^ ( 77 , T, N, M, 0,0) 


1 


(2tt) 


'[0,2ir)x[0,27r) 


dfadfa 0-^2 ,N \fa 


< 


(2tt) 


'[0,27r)x[0,27r) 


d(p\d(t>2 


d{^ N \rj,T,^,M,0,0) - 


Therefore, in order to prove Theorem [TJ it is sufficient to show that the 
integrand in the right-hand side is vanishing in the multiple limit in (14.121) . 
The integrand is estimated as 

d£’ N \r,,T,X,M,0,0) - a£ N \cj>) 

< d£' N) ( 77 , T, N, M, 0,0) - a$ N) ( 77 , T, N, M, fa, fa) 

+ d[2 ,N \v,T,J^,M,fa,fa) - d[2’ N) (0,oo,'N,M,fa,fa) 

+ <5-i2^ ,7V) (0, 00 ,3N, AT, ,0 2 ) - a[ 2 ,N \<j>) . 


As shown in (110.11) above in the present section, the first term in the right- 
hand side is vanishing in the infinite-volume limit A /* Z 2 . Relying on the 
estimate m, we can show that the second term is vanishing in the double 
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limit r) —>• 0 and T —> oo after taking the infinite-volume limit. Finally, the 
result (18.21) implies that the third term in the right-hand side is vanishing 
in the limit. 

Since the above Hall conductance d[ A,N ^ ( 17 , T, N, M, 0,0) with the vanish¬ 
ing phases is equal to the Hall conductance ,N ^ ( 17 , T, (N, M) in the case 
of the periodic boundary condition in the right-hand side of (14.12p . we have 
proved the fractional quantization (12.21) for the Hall conductance a \2 in the 
limit. 


Appendix A. The Correction to the Hall Conductance 
In this appendix, we show that the correction 

4or ■■= J ^ds i]se vs u^’ N \[x {A \T M (k,£)- s), J^ } (k,£)]) 

in (|4.10|) to the Hall conductance is vanishing in the limit r] —>• 0 after taking 
the limit T —» 00 . 

We denote the ground-state vectors of the unperturbed Hamiltonian 
by $ 0 ^ with the energy eigenvalue m = 1 , 2 ,... ,q, and the 

excited-state vectors by <3?^ with the energy eigenvalues E^ N \ n = 1 , 2 ,.... 
Using the definition (I4.8j) of s), one has 


4 A,JV) ([x (A) (r M (M); s), (M)]) 


m= 1 n> 1 L 


— C.C. 


Further, one has 


(A ,) /;, 


, ris ifE^-E^hs 

as se 1 e y °’ rn ' 


iT 


e {N) _ (N) , 

-^0 ,m + L '\ 


e ~vT e i{E^-E^)T + 


[ h n ~ E Q m + IT] j 
_ e -vT e i{EL N) -E< N 2 )T' 


[ h n - h 0 ,m + l d) 

From these observations, we can evaluate the integral as 

|x (A) (r M (M))|| || 4 1 ) (M)||, 


lim \I COT \ < -^ 2 |I^ (A) ( 
T— >00 AE 2 


where we have used the Schwarz inequality and the assumption of the spec¬ 
tral gap, E^ — > A E. Thus, we can obtain the desired result in the 

limit 77 —>• 0 for M and IN’ which are finite. 
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Appendix B. Proof of Theorem [2] 

B.l. Constructing the low-energy excitation. Let us consider a generic 
Hamiltonian H^ with finite range hopping amplitudes and finite range in¬ 
teractions on the lattice A. We assume that the Hamiltonian H^ commutes 
with the total number operator of the fermions. Let N be the number of 
fermions, and fix the filling factor v = A/|A|. We denote by H AA) the 
restriction of H A) onto the eigenspace of the total number operator with 
the eigenvalue N. Further we assume that there exists a uniform spectral 
gap A E above the lowest q eigenenergies Em = 1,2 ,... ,q. We write 
for the corresponding q eigenvectors, and define 


SE := max 

m,m' 


(I 


E, 


(N) 
0 ,m 



We do not assume that the q energies E ( ^J are (quasi)degenerate. Therefore, 

we allow SE > 0 uniformly in the size |A| of the lattice. We write e[ N * for 
the q + 1-th eigenenergy of the Hamiltonian H AA) f r0 m the bottom of the 
spectrum. Let be the projection onto the sector which is spanned by 

the q vectors and define the expectation for the sector as 


u 


(A ,N) 
0 


(...) 


1 

q 




(N) 

0,m’ 


m= 1 



Proposition 6. Suppose that, for any given small £ > 0, there exist a local 
observable ao with a compact support, and a positive constant co such that 

(B.l) (<hg ) ,aSP £/2 a 0 $g ) )>co>0 

for a vector in the set of the q vectors {Tq^}, where P e / 2 is the pro¬ 
jection onto the energy interval [e[ N \ e[ N%> +e/2]. Here both the support 
size of the observable ao and the constant c can be chosen to be independent 
of the size |A| of the lattice. Then, there exist a local observable a with a 
compact support and a positive constant Vq such that 


(B.2) 


A E < 


OJ 


( A A0( a *(i _ P ( A ’ N ^[ H W,a}) 

4 A ’"Va-Po (A ’"V) 


< A E + 2 SE + 2e 


for any A satisfying |A| > Vo. Here the size of the support of the observable 
a is independent of the size |A| of the lattice. 


Remark . If the condition (IB.ID does not hold, then there must exist another 
infinite-volume ground state which cannot be derived from the q vectors. In 
other words, all of the low energy excited states are derived from a local 
perturbation for the ground state. 


In the rest of this subsection, we will prove the upper bound of (IB.2D 
in Proposition [6] because the lower bound can be easily obtained from the 
definitions. 
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We decompose the energy interval in Proposition [6] into M small intervals 


as 


M 


[E[ N \E[ N) +e/2] = U {E{ N) + (n-l)s/(2M),E[ N] + ne/(2M)\, 


n =1 


where M is a large positive integer. Then, there exists a small interval 
[E[ N) + (k - l)e/(2. M),E[ N) + ks/(2M)] 

such that 


(B.3) 

where k is an integer in {1,2 ,,M}, and P*j 2 is the spectral projection 
onto the small interval. 

Let K be a large positive number satisfying e\[K f 4 < M < ey/K/4 — 1, 
and write 




co 


and 


A E := E{ N) - E^O, 


R, := {k- 1 / 2 )-—. 
y ' ’ 2M 


We use the idea of an energy filter [6]. We define three operators as follows: 


a o 


1 oo 

L= / dt e itHW a 0 e~ itHW 

!t tK y_ oo 


and 


5o(^i) ; — 


a(T x ) := 


\J2 txK J— 
i r +T i 


\Z 2 ttK J—7 


dt e^ooe^e^+^e-^, 


f+Ti 


\/2ttK J 


dt e itHW a 0 e~ itHm e -i^+k)t e -P/( 2 K), 


Ti 


where the cutoff Xi is a large positive number, and is the restriction of 
the Hamiltonian H to the region UcA, We choose H so that supp ao C 
fl. Clearly, the operator o(Ti) is local and has a compact support in H. We 
write a = a(T\) for short, and we will prove that this observable a satisfies 
the upper bound of (IB.21) in Proposition [6] for appropriately choosing the 
parameters. 

We decompose the projection (1 — Pg A,JV ' ) ) into two parts as 


1 - P, 


(A,IV) _ 


0 


Plow T Phiglu 


where Pi ow is the spectral projection onto the interval [pj^, e[ N ^ + e + 5E\, 
and Phigh the spectral projection onto (e[ N * + e + SE, +oo). Then, we have 

(B.4) 4 A,A °( a *(l - T > Q A ’ Ar) )[E^ A \ a]) 

= 4 A,7V) («^ l0 w[Tr (A) ,a]) +4 A,iY) (a*T> high [iL (A ),a]). 
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Since supp a C fl, we can find a region SfDllso that [H^ A \a] = [i4^\a]. 
Using this and the Schwarz inequality, the second term in the right-hand 
side of (IB. 41) is evaluated as 

(B.5) 

4 A,JV) (a*P h igh[^ (A) ,a]) = |4 A,iV) (a*Phigh[^ (f2) ,a])| 

< \Jub A ' N) (a*fhi g ha)wo A,Ar) (> a l* > a ]) 

< 2\\H^\\ \\a\\^4 A ’ N \a*P high a). 

From the definition of a = a(Ti), we have 

i r+T\ 

||a|| < —j= / dt ||a.o||e — * /( 2jF4r ) < ||a 0 ||. 

_ VzkK J-Ti 

Substituting this into (IB. 51) . one has 

(B.6) co^ N \a*P high [H^\a]) < 2\\H^\\ ||a 0 || 

The first term in the right-hand side of (IB. 41) is written as 
4 A, " ) (n*fl„ w [.£r< A >,a]) 

= \ E(fS.o*flo,[fl (A) ,o]fS> 

y m=l 

= - a t £ <<4 «**<''>> <*'«.tg)(4*>-C). 

m ~ l E^ N) e[E ( 1 N) ,E ( 1 N) +£+5E] 

where we have written for the excited-state vectors with the eigenen- 
ergies E^\ Therefore, we obtain 

(B.7) J 0 A ’ N \a*P low [H^ A \a]) < (AE + 25E + e)J A ’ N) (a* P low a). 

Clearly, the denominator in the expression ([B.2D of the excitation energy 
can be decomposed into two parts as 

4 A,A V(1 - = ^V^iowa) +4 A ’ Ar) (a*P high a). 

From these observations, it is sufficient to evaluate the first and the second 
terms in this right-hand side. 

First, we will show that the operator a is approximated by the operator 
do- Since one has 

a — do = a — do(7i) + do(7i) — do, 
we will show this right-hand side becomes small. Note that 

d 0 - do(Ti) = -pL= dt e itHW^ e -itHW e -i(AE+!i)t e -^/(2K) 

v27Tj K JTi 

+ -=L= [~ Tl dt P tHW a 0 e~ itHW 
V2t tK J -oo 
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Therefore, we have 

(B.8) ||o 0 - a 0 (Ti)|| < -^j== J e~ t2/{ - 2K ^ < ||a 0 ||e _T i /(2X) . 


We write 


and 


a[ K \t) := e itH(A) a 0 e~ itH(A) 
4 n) (t) :=e itH(Q> a 0 e- itH(n) . 


We have 

5o(Ti) - a = —L= j + J‘ dt l4 A \t) - 4 n> (t) ]e -i(A«+s)< e -«V( 2 K) 
The norm is estimated as 

1 /»+Ti 

(B.9) ||ao(Ti) - a|| < j ^ dt e~ t2 ' {2K) \\a ( 0 A \t) - a^ 0) (i)||. 

The norm in the integrand is estimated as (22j 

(B.10) ||a[, A) (t)-a[> n) (t)|| < sgnt [ ds ||[(iL (A) - H^),a^\t - s)]||, 

■Jo 


where 


sgn t := 


< o, 
w 1 ’ 


t > 0; 
t = 0; 
t < 0. 


The proof of the above inequality (IB.lOp is given in Appendix [Cj 

In terms of local operators hx with the support X, the difference between 
the two Hamiltonians are written as 


tf(A) - H (n) = Y h * 

x 

xn(A\Q)^0 


Therefore, we have 

(B.ll) \\[(H^-H^),a^\t-s)}\\< Y ||[^,4 n) (t- S )]||. 

a: 

xn(A\Q)nn^0 


In order to estimate the summand in the right-hand side, we recall the Lieb- 
Robinson bounds EUB El E3J: Let Ay,Bz be a pair of observables with 
the compact supports, Y, Z, respectively. Then, the following bound is valid 
0: 

(B.12) || [Ay(t), Bz] || < C|| Ay || ||^|||y||Z| exp[—/x dist(Y, Z)\ [e v ^ — 1], 

where Ay(t) = e ltH<X ' Aye^ ltH< '' > for the Hamiltonian H^ with finite-range 
interactions, and the positive constants, C, v and /jl, depend only on the 
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interactions of the Hamiltonian and the metric of the lattice. Using the 
Lieb-Robinson bounds, the right-hand side of (IB. Ill) can be evaluated as 

< Const. ^2 |suppao| ||ao|| (e v ^~ s ^ — l)e -/ir 

a: 

xn(A\n)nJV0 

< Const.|supp ao| HaolKe^* - ^ — 1 )Re~ flR , 


where r = dist(X, supp oo), R = dist(supp ( H — H^), supp ao), and the 
positive constants, fi and v, depend only on the parameters of the present 
model. Combining this, (IB. 911 and (IB.10D . we have 

||ho(Ri) — a|| < Const.|supp ao| ||ao||Rexp[nTi — /j,R]. 

Combining this with the inequality (IB. 81) . we obtain 

(B.13) ||a - a 0 || < ||a - a 0 (Ti)|| + ||a 0 (Ti) - a 0 || 

< Const.|supp ao| ||ao||Rexp[wTi — fiR\ + ||ao||e _T2// ( 2A ^. 


Thus, one can make the difference between a and ao small by appropriately 
choosing the parameters, T\, R and K. 

Next, we prepare some estimates about the operator ao- We denote by 
Qri' 1 ' 1 the excited-state vector with the energy eigenvalue E^ > e[ N%> for 
n > 1. The matrix elements of ao between the excited states and the lowest 
energy q states are computed as 

1 Z* - |“00 

= (^\a 0 ^)e W [-K(^ m ) 2 /2], 
where we have written 


^n,m 


:=EW-EW-AE-k. 


Using this relations, one has 

(•ffl.a5Aw.aoO 




K \2 
n.m) 


E^ n) >e[ n) +e+SE 


< E 

E^ N) >e[ N) +e+5E 

— /a,W) _* p _ Ke 2 /4 

— \^0,m! a O'niighao < f > o,m/ e 


Immediately, 

(B.14) ^ A,Af) (aoRhi g ha 0 ) < ^ A,JV) (aoPhi g hao)e _A£2/4 < ||a 0 || 2 e _A£2/4 
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Similarly, we have 
q X ^ A,Ar) (aoPiowOo) 

= ^(C^owaoO 


m= 1 
Q 


E E 

m ~ 1 e^ n) g[e[ n \e[ n) +£+5e] 
<? 

E E 

m ~ 1 Ei N) e[E[ N \E[ N) +£+8E] 


<<2,«sei<euc>e- A<£ - )a 


e N) ,e[ N) +e+8E] 


For 


EW € [£f° + (k - l)e/(2M),£f° + «e/(2M)], 


we have 


. < JL 

4M “ “ AM ' 


Combining these observations, the inequality (IB. 311 and the assumption M > 
sVK/A, we obtain 


(B.15) 


■PlW.) a W.l'.oSP^oofiy)— > 




(iv) \ ^ 1 ^ e 1 c 0 


Ti = —=Ke, and riR = —(2\/2w + e)Ks. 
y/2 4 v 


We choose the parameters, Ti, JC and T so that they satisfy 
(B.16) 

Then, the inequalities, (IB.1311 and (IB.14H . are, respectively, written 
||a — ao|| < Const.|supp ao| \\ao\\Re~^ eR 

and 

w o A ’ JV) (“o- p higha 0 ) < ||ao|| 2 e _M£jR , 
where we have written 


H := 




2y/2v + £ 

Using these inequalities, we have 
(B.W) 

uj^ A,N \a*Ptigta) = ujQ A,N \(a - a 0 )*P hig ha ) + ojl A,N \a* 0 P high (a - a 0 )) 
+ w^ A,Af) (aoPhighOo) 

< 2||a 0 ||||a - a 0 || + ^’^(agThighao) 

< Const.|supp ao| Haol^Te - ^^. 
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Similarly, 

(B.18) \4 A ’ N \a*P low a) - 4 A ’ N) (a 0 P low a 0 )\ 

< \uj { 0 A ' N \(a - a 0 )*Pi ow a)| + |w^ A,JV) (a 0 Piow(a - a 0 ))| 

< 2||a 0 ||||a - ao|| 

< Const.|supp ao| ||ao|| * 2 i?e _ ^ £R . 

Now let us estimate the excitation energy. Using the inequalities, (IB. 61) 
and (IB. 71) . we obtain 

4 A,JV) (a*(! - P 0 (A,JV) )[i7 (A) ,a]) 

4^(1-P^)a) 

_ 4 A ’ N) (a*P low [H^\a]) +4 A ’ N \a*P high [H^\a}) 

4 A,N \a* Plow®) + 4 A,N \ a * P high a) 

„ (A E + 2 5E + £ )4 A,JV) (a*^owQ) + 2||i7^)|| ||a 0 || 

4 A ’ N \a* Piowa) + 4 A ’ N \ a *Pi iigh a ) 

We can choose the region (l so that 

\\H^\\ < Const.i? 2 4 

for a large R. Combining these observations with the inequalities ()B. 171) 
and (IB.181) . we obtain 

miql 4 A,jV) («*(i-^o (A,jV) )[g (A) ,«]) 

‘ ^(a-Cl-P^Ja) 

< (A E + 2 5E + £)J A ' N) {a* 0 P low d 0 ) + 0(R 5 / 2 e~^ R / 2 ) 

4 A ' N) («Sflowd 0 ) + 0(Re~^ R ) 

for a large R. In order to estimate the right-hand side, we write the inequal¬ 
ity (IB. 151) as 

J A ’ N) (a* 0 P low a 0 )M > > 0. 

Combining this bound, the assumption M < ey/~K/4 — 1 and (IB. 161) . the 
right-hand side of (IB.191) with a large R yields the desired upper bound (IB. 2D 
for the excitation energy. 

B.2. Twisted phase dependence of the excitation energy. We denote 
by Hq A \cJ)) the unperturbed Hamiltonian of the present system with the 
twisted boundary condition in the first direction with the angle 4>- We write 
H^ A ' N \cj)) for the Hamiltonian restricted onto the sector of N fermions. 

We denote by PQ A,N \(j)) the spectral projection onto the sector which are 
spanned by the q vectors with the lowest q eigenenergies. We denote by 
^*0m(^) the low-energy vector with the eigenenergy m = 1,2,... ,q. 

The expectation value in the low-energy sector is given by 
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Although we have assumed that the ground state for <f> = 0 is g-fold degen¬ 
erate, we do not necessarily assume that the low-energy sector for (f) ^ 0 is 
degenerate. 

Let a be a local observable with a compact support. The excitation energy 
due to the local perturbation a is given by 

aTVa - p t) A ' N) A) (0)>a];<ft) 

(«*(! - p o (A,iV) (<£))«; 4 

We want to show that this excitation energy is almost independent of the 
phase (j) for a large lattice A under the assumption of the nonvanishing 
spectral gap above the low-energy sector. 

The numerator of the excitation energy is written 


4 A,A ° i a * (! - p o A ' N) (<£)) i H o A) 0 )>«]; 4>) 

= ( a * i H o A) (</>), a]; </>) - 4 AlJV) (a*P 0 (A,JV) O) [H { 0 A) ((f)), a]; 0) 

The second term in the right-hand side is written 

4 A,A ° (a*-P 0 (A,Ar) (0) [ F o A) 0£)> °]; 4 >) 

= 1 -it [4 A) ,°]<2m> 

’ m=l m'=l 


m,m' 


in terms of the low-energy vectors 4/m ( ( f > ). Therefore, we have 
(B.20) 

|4 A ’ W) (a*Po A,N) (</>)[H { 0 A) ((f)), a];4>)\ < 8E(4>) u^ ,N) (a*P^ A ’ N) ((f>)a; <f >), 
where 

5E((j)) := max{| E^(<f>) - E^ty) |}. 

m,m' 

Consider the excited states with the energies which are larger than the 
maximum value, maxi<j< (? {£ , Q4 (</>)}, of the lowest q eigenenergies. We 
denote by <&n*\<f)) the excitation vector with the eigenenergy En N \cf>) for 
n > 1. We take the subscript n of the eigenenergies {E^\(f>)} to satisfy 
E[ N \(f>) < E^\(f)) < .... We assume that there exists cj>o € [0,27r] such 
that 
(B.21) 

min {e\ N \(/))} — max max{£^ (<f )),... ,E^ N \(j))} > AE mm > 0, 

4>e[o,M </>€[ o,0o] u ’ u ’ 9 


where AE mm is independent of the size |A| of the lattice. That is to say, 
there is a uniform lower bound for the spectral gap above the low-energy 
sector. 

The denominator of the excitation energy is written 


OJ, 


(A,TV) 


(a*(l- PQ A ' N \(j))) a ;<t>) = J 0 A,N) (a*a; <f>) - u)^ A,N) (c 


*P< 


(A,TV) 


0)a;</>). 
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The second term in the right-hand side is written 

(B.22) ^ A ’ N) (a*Pi A ’ N) (^a^) = -TV a* P^ N \^)aP^ N) ((/>). 


Under the above assumption on the spectral gap, the projection P^ A,N \(j)) 
can be written 


P 


(A,TV) 


((j>) = ——: ® dz 


2n i 


z — H, 


i A) « 


on the sector of N fermions for cj) € (0, 0o] - This is also differentiable with 
respect to 0 as 


d_ 

d(f> 


i KN) w = 


1 

2iri 


dz- 


1 


z-H {a \4>) L 

As we showed in Sec. 0 the operator 


d(\) 


1 


z-4 A) W 


B(4>) := ^ 


(A) 




has a support which is localized near the boundary. By differentiating and 
integrating (IB.22[) . one has 

4 A,JV) ( a *^o A Ar) ( l ^o)a; 4>o) ~ 4 A,JV) (a*-P 0 (A,Ar) (0)a; 0) 


- / 0O #4 

q Jo # 


TVo*P„ (A,JV) (^)aP ( (A,JV) 


o o (0)- 


We will show that the integrand in the right-hand side is small for a large 
volume. In consequence, the difference of the expectation values in the left- 
hand side is almost independent of (f>o for a large volume. 

Relying on the differentiability of the projection operator, one has 

(B.23) -^-Tr a*P^ A ’ N \(j))aPQ A ’ N \(j)) 


= Tr a* 


d 


aPQ A,N \cj)) + Tr 


= Tr aP^^—P^^) + TV a*P o (A>JV) (0)a— P^ N \cf>). 

Since the first term in the right-hand side is obtained by interchanging a* 
and a in the second term, we will treat the second term only. The second 
term is written 


Tr a* P, 


(A,TV) 


d 


o 


W “^ P o 


(A,TV) 


( 0 ) 


1 

2tti 


dz Tr a* P 


(A,TV) 


0 


(</>)«- 


z-H, 


i A, (0) 


-*>«)- 


1 


z-H, 

i(JV) 


i A V) 


- 7Tr\ - 7Tr\ - +C.C. 


m= 1 n> 1 

Using the integral identity 




r c 
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the first sum in the right-hand side is written 

poo 9 _ 

ds T, 

■' 0 m=ln>l 

x exp [~{E^\(f>) - ^2(0))s] 

per Q _ 

"'° m=lra>l 

x exp[-(E^(^>) - E^((j)))s} 

poo Q _ 

- / ^ ]T ^<^(^),a*P 0 (A ^ ) (^)a^(^))($W(^),5^)^(^)> 

•J cr „—i 


x exp[-(.E^ 0) - #o2(0)) s ]» 

where c is a positive constant, and r = dist(supp a, supp B (</>)). The second 
integral in the right-hand side is estimated as 


poo 9 _ 

/ ds J2 £<O*)’ a * P 0i A,JV) («fl$i JV) («){$i JV) (« I B(«CW) 

Jcr __i „ 


m=l n>l 


x exp[-(£;W(^) - £$J2(0 ))s] 


< 


AU E 


|a|| 2 115(0)|| exp[-A£ min cr] 


where we have used the assumption (IB. 2111 on the spectral gap. Since 
||i?(</>)|| = 0(L( 2 )) and r = dist(supp a, supp B{(f))) = 0(LW), the corre¬ 
sponding contribution is exponentially small in the linear size of the system. 
The first integral in the right-hand side is written 


(B.24) 

where 

with 


per 

~q ds u}^' N \a*P^ A ' N \(j))aB l ' A ' > (is; <t>)\4>), 
Jo 


B^\z- (/>) := for z€C 


(A), 


B{^) := B(</>) - P^ N \(t>)B(ct>)P^ N) ^). 

Since r = ©(L^), it is sufficient to show that the integrand of (IB.2411 is 
exponentially small in the size L^ l \ For the integrand of (IB.2411 . we write 

f(z) = ujQ A,N \a*PQ A ’ N \(j))aB^ A \z] (/))-, (/)) for z € C. 

Using the contour integral in the complex plane, one has 


(B.25) ul A ' N \a*P^ N \(j))aB {A \is-,(j))-,(j)) = £> ^ dz 

2m J z — is 


1 

2ni 


rT ' 2 f(t) , 1 

J — dt + 


m 


'-t 2 


t — is 


2ni 


■dz 


z — is 


for s > 0. Here, the first term in the right-hand side of the second equality is 
the integral along the real axis from — T 2 to T 2 with a large positive number 
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T 2 , and the second term is the integral along the semi-circle z = T 2 e ld for 
0 E [0,7r]. The second integral is estimated as 

- 2 AP min r 2v /i - 2 s/T 2 

for 2s < T 2 . The proof is given in Appendix [Dj We recall that ||P(</>)|| = 
O(L^) and that 0 < s < cr with r = O(L^). Therefore, the corresponding 
contribution is vanishing in the limit T 2 —» 00 . 

In order to evaluate the first integral in the right-hand side of the second 
equality of (IB.251) . we use the technique in 0 01231. Following them, we 
write 

f(t)e as2 = f(t)e~ at2 + f(t)(e as2 - e~ at2 ). 

By definition, 

/(f) = uj { 0 A ’ N \a*P^ ,N \ct>)aB {A) (t- <f>)\4>). 

Using the identity, 

aB (f; <f>) = B^ (t; 4>)a + [a, B (A ) (t; <fi )\, 

one has 

fit) = 4 A,A ° ( a *-P<S A,N) (</>)B (A) (t; <f>y, 0) 

+ u [ 0 A ' N \a*P^ A ’ N) ((l))la, P (A) (t; </>)]; <j>) 

= w o A,7V) (^ (A) (*; ^ ap r ] (</>)«*; 0) 

+ ^(a*P 0 ^ W [a, -B (A) (t; </)]; </) 

- 4 AlJV) (a*T > Q A,Ar) (</>) [a, P< A ’"> (/)f?( A ) (f; (<«]5 <^>), 

where 

B^\t-(f>) := e iH ° A)t Bi(j))e- iH o A)t 

for t E R. Consider the third term in the right-hand side of the second 
equality. Except for the factor q _1 , the term is written 

Tr a*P^ A ’ N) (ft [a, P^’ N) (^B^ (f; c/>)P^ N) (</)] P 0 (AliV) (cj>) 

= TV a*P 0 (A ’ JV) (</)aP 0 (A ’ iV) (</)p( A )(f; cj>)p( A ’ N) (d>) 

- Tr a*P 0 (A ’ JV) («/)p( A )(f; $P 0 (AlJV) (cj>)ap( A ’ N) (0) 

= Tr a*P 0 (A ’ JV) (</)aP 0 (A ’ iV) (</)P (A) (f; $P 0 (A>JV) ($ 

- TV ^ A ’ N) (0). 

Clearly, by interchanging a* and a, this contribution changes the sign. 
Therefore, the contribution is canceled by the corresponding contribution 
coming from the first term in the right-hand side of the second equality of 
(IB.231) because the first term is obtained by interchanging a* and a in the 
second term in (IB. 231) . 

From these observations, in order to evaluate the first integral in the 
right-hand side of the second equality of (IB.231) . it is sufficient to estimate 


(B.26) 


1 

2m 


m 


■dz 


z — IS 
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the following three integrals: 

(B.27) h := F dtcj>)ap( A ’ N) (0)a*; </>)e 
27T* J_t 2 t — IS 

(B.28) / 2 := -L F dt 4 A,iV) (o*P 0 (A,Ar) (^)[a, # (A) (t; $]; ^e -0 * 3 

Z7TZ J_ r j 1 ^ Z — 28 


(B.29) J 3 := ^ F dt l^L{e as2 - e~ at2 ) 

2m J_t 2 t — is 

The first and second integrals can be estimated as follows: 

(B.30) |/i| < l|a||2|l ^ WI1 exp[-(A£ min ) 2 /(4a)] 

in the limit T 2 —>• 00 , and 

(B.31) |I 2 | < \\a\\ 2 \\B{(t))\\( ~^=\e~ aS2 + Const. |supp o|e ,,s - /ir 

\ y/'KOt O 

Here, S' is a positive number; v and /r are positive constants which are 
determined by the model’s parameters. The proofs of (IB.301) and ()B.31D are 
given in Appendices [E] and [Fj respectively. As to the third integral with 
0 < s < cr, we choose a = AE mm /(2cr). Then, one has 

(B.32) |I 3 | < i||a|| 2 ||B(^)|| e - A£min W2 

in the limit To —>• 00 . The proof is given in Appendix [Gj 

For the same a = AE mm /(2cr), the upper bound (IB.301) for |/i| becomes 
the same as that for | T 31 as 

IAI < \\\af\\B(4,)W^ E ‘"^. 

As to the upper bound (IB. 311) for |1, we choose S = cr and 


v + AE min /2 

with the same a = AE mm / (2 cr) as in the above. Then, we have 


\h\ < ||a|| 2 ||i?(</>)|| ( —== + Const.|suppa| )e ^ r , 
yjixfir 


where 


A E D 


V = 


2v + A E' 


— u. 

Tin • 


Since AE mm c/2 = fi, we have 


\h\ + \h\<\\a\\ 2 \\Bm^ r . 

We recall ||H((/>)|| = O(L^) and r = dist(supp a, supp B{4>)) = ©(l/ 1 )). 
Therefore, all of the above integrals are exponentially small in the linear 
size ZA 1 ) of the system. As a result, we obtain 


(B.33) |4 A,A ° {a*P^ N) Oo)a; fo) - J A ’ N) {a*P^’ N) (0)a; 0) 

< Const.|supp a|||a|| 2 ||B(^>)||e _ ^ r . 
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Clearly, we can write 

(B.34) 4 A ’ iY) («^o A ’ 7 V) (^o)a;</>o) = u {A ’ N) (a* P^’ N) (0)a; 0) + ei(A), 

where ei(A) is the exponentially small correction. 

Since only the hopping terms near the boundaries in the Hamiltonian 
H^ A \(j)) depend on the phase cf>, the commutator (</>), a] is independent 
of the phase (j) and becomes a local operator with a compact support which 
is independent of the lattice A. Therefore, we obtain 

(B.35) u { 0 A ’ N \a*[H { ^\(j) 0 ),a}-,(j) 0 ) = u^ A ’ N) (a* [H { Q A) (0), a]; 0) +e 2 (A) 

with the small correction e 2 (A) in the same way. Furthermore, we have 

(B.36) Uq A ' N ^ (a*a~, <j> 0 ) = u}q A,N ^ (a*a; 0) + e 3 (A) 

with the small correction 63 (A). 

We define 


SE max := max \8E((j))}, 

4>e[o,4>o] 


and 

A E(<t>) := E[ N \(j)) - rnax{^ } (</>),... , E$ (0)}. 
We assume that there exists a small e > 0 such that 


(B.37) A E min - 8E max > e. 

Let us consider the Hamiltonian Hq A \(/)q) with a fixed 4> o- We require 
the same assumption as in Proposition [ 6 l For the same e as in the above, 
we can find a local operator a which satisfies the excitation energy bound 
(B.38) 


a4 A,Ar) (a*(l - P^ A,N \^ Q ))[H^ A \^ Q ),a];^ Q ) 


< AE((f>o) + 25E{<j)Q) + 2 e 


oj { 0 A,N) (a*(l - P 0 (A ’ Ar) ( <? i 0 ))a; <fo) 
as in Proposition [ 6 l We write 

6 4 (P) := J A ’ N \a*P^ N \cl )0 )[H {A \ ( l )0 ),a}-A^ 

for short. The same argument in the subsection IB. 1 1 yields 

(B.39) 4 A,JV) («^o A,iV) (^o)a;</>o) = 0(Re~^ R ). 

Combining this with the bound (IB.201) . we have 


e 4 (P)| < SE(</>o)ui A ’^{a*P^{<h)a\<h) = 0(Re~^ R ). 
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From (IB.341) . (IB.351) and (IB. 361) . the excitation energy for the Hamiltonian 

Hq A \4>o) is written 

(B.40) 


u q A, n ' > ( q *(i ~ 4 A,A 4^o))[#<4(^o)> Q ];^ 

w 0 A,iV) ( a *( 1 - p o (A,iV) (^o))a;^o) 

_ Vq A,JV) (q* [ h q A) Oo), a]; 4>o) - Uq A ’ N) (a*PQ A ' N) (4) [Hq A) Qo), a]; 0o) 

4 ^ (°*a; O) - 4 A,JV) {a*Po A,N) (0o)a; 0o) 

4 A ’ JV) (a* [ff< A) (0), a]; 0) + 6 2 (A) - e 4 (fl) 

4 A>JV) (a*a;0) +e 3 (A) - u;^’^ (a*P 0 (A,iV) (0)a; 0) -ei(A) 

We assume that the ground state of the Hamiltonian Hq A,N \ 0) with 0 = 0 
is (quasi)degenerate, i.e., 


SE( 0) -A 0 as |A| —»• oo. 

Under this assumption, from (IB.201) . we have 

WQ A,Ar) (a*PQ A,iV) (0)[HQ A) (0), a]; 0) -A 0 as |A| -A- oo. 

Further, from the same argument in the subsection lB.il (IB.341) . (IB. 361) and 
(IB.39I) . one can show that there exist a sufficiently large |A|, a sufficiently 
large R and a positive constant c' 0 such that 

J A ’ N \a*(l-P^ N \o))a-,0)>^, 

where the positive number M is defined in the subsection IB. II From these 
observations, the excitation energy of ()B.40D can be estimated as 

4 A,iV) (a*(l ~ ^o (A,7V) (0o))[^p A) (0o) ? q];0 o) 

4 A,A) («*(! “ 4 A,iV) (^o))a;0o) 

4 A,iV) («*(l - P^ N) m[H {A \0),ay,0) 

4 A ’ N) {a*(l-P^ N \0))a-0) 

> AE(0) - e 


for a sufficiently large |A| and a sufficiently large R. Combining this with 
the upper bound ()B.38D of the excitation energy, we obtain 

(B.41) A£(0 o ) > AE'(O) - 26E(4>o) - 3s. 

Thus, if AE(0) > 25E((f>o) in the infinite-volume limit, then the excitation 
energy for 0o is strictly positive, i.e., AE((j>o) > 0. Further, if the splitting 
SE (0o) of the energies in the low-energy sector is vanishing in the infinite- 
volume limit, then we have the bound Ai£(0o) > AE(0). 


B.3. Twisted phase dependence of the averaged energy. Consider 
the same situation as in the preceding Section IB.2I But, instead of the 
energy gap condition ()B.21D . we consider more generic condition, 

(B.42) 

min {E[ N \(j))} — max max{£'o^(0),... ,Eq N \(/))} > AE mm > 0, 
1,02] 0 e[</>i,</> 2 ] ’ ’ 
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for the interval of the phase <j>. Here, we assume that the positive 

constant AE min is independent of the size |A| of the lattice. 

We decompose the Hamiltonian (cp) into two parts as 

^ A) (0)=^ A+) (0) + ^ A " ) (^) J 

where the supports of the operators H {] ; (</>) are included in the regions, 

A + = (j-L (1) /2,-L (1) /2 + «o] U [—i? 0 ,i (1) /2]) x [-L (2) /2, L (2) /2] 

and 

A' = ([~lW/2,R 0 \ U [-fio+f (1 ) /2,f (1) /2]) x [-LW/2,lW/2], 

respectively. Here, i?o is a positive constant which is independent of the size 
|A| of the lattice. Clearly, we have 

(B.43) 4 A ’ JV) (F ( S A) (0);^) =4 A ’ Af) (^ A+) (0);^)+4 A ’ JV) (^ A " ) (^);^). 

Consider the first term in the right-hand side. As shown in Sec. [21 one can 
change the position of the twisted hopping amplitudes of the local Hamilto¬ 
nians in the total Hamiltonian by using the unitary transformation. There¬ 
fore, we can find the unitary operator U- such that in the transformed 
Hamiltonian 

the local Hamiltonians with the twisted hopping amplitudes are arranged 
along the center line with the first coordinate ~ —L ^/4 in the region 
A - . We write 

^ A+) := U^H^ A+ \(j))U-. 

By definition, this Hamiltonian Hq A _ ^ is independent of the phase </>. Fur¬ 
ther, one has 


(B.44) dist(supp H^ A _ \supp B_(0)) = 0(L (1 ^), 

where the operator 

comes from the twisted hopping terms in the Hamiltonian. In other words, 
we choose the unitary transformation U— so that the condition (IB.44[) is 
satisfied. 

The first term in the right-hand side of (IB.4311 is written 


Wn 




1 . 


(A+) „(A ,N) 


=-TiH^_>P 0 


1 1 

q 2ni 


dz Tr H. 


(*) 

(A+) 


0 ,- 


z — Hq A } (4>) 
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Relying on the assumption (IB.42D on the spectral gap, one has 


^4 A ’ 7 V) (^o A+) ( < /’);^) 

1 1 

q 27 vi 


dz TV //, 


(A+) 


°’ z-H^) 




z-H^) 


By integrating the left-hand side, one obtains 


w (A,A 0 (^ A +)(^ 2 );^ 2 ) - w J A,JV) (fl'i A+) (^ 1 ); 0 1 ) 

= C «^ A,A ' ) (^ A+) W!'#')- 

Combining these with (IB. 441) . in the same way as in the preceding Sec¬ 
tion IB.21 we can prove that this left-hand side is exponentially small in the 
linear size L ^ of the system. Since the second term in the right-hand side 
of (IB.431) can be handled in the same way, we can prove that the difference 

4 a,jv) (^ a) (^ 2 ); ) - w£ a,jv) (-Ho A) 0 £i); ) 

is exponentially small in the linear size L ^ of the system. 

Next consider the situation that the degeneracy of the energies in the 
low-energy sector is lifted. Namely, we consider the situation that there is 
a spectral gap between their energies in addition to the above assumption 
(IB. 4211 on the spectral gap. 

Let mi,m 2 be integers satisfying 1 < mi < m 2 < q. We assume the 
existence of the nonvanishing spectral gaps as 


min rnin{P ( 


0,7711+1(0)’ 


’ -^O ,7712 


m 

(N), 


— max max{Lj x ( 0 ), 
</>e[ 0 i, 02 ] 


A N 2,m ^ AE - in >° 


and 


min min{45 2 + 1 (0), ■ ■ ■ > 4^(0)} 

(p£[(pl,(t>2\ 

~ . max max{42 1+ i(0)’ • • •, 42, (0)} > A ^+ m > 

0 G[ 01 , 02 j 

where A£“ m and A E™ m are independent of the size |A| of the lattice. 
When m 2 = q, the second condition is replaced with the condition (IB. 4211 . 
Let P ( j A y !V ) (0) be the projection onto the sector spanned by the eigenvectors 
with the eigenenergies (£q^(0 ),... ,£^^(0)}, and P^ N \(j)) the projec¬ 
tion onto that with (£q^ i+1 (0), ... ,Pg^] 2 (0)}. We write 

4 A 1 N) (--- ; 0) = —Tr (• • • )P O (A 1 JV) (0), 

mi 

and 

4 A w’(■" ;« = ("') K+^W) + dPV)] ■ 
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Then, one has 


, X A ’ N )ru 
w o,+ l-'u 


o A \fafa) 

m 2 


m 2 — mn\ 


(A ,N) / rr(A) / i \ . I \ _ 

w o,tot (^o (fa fa) 


mi 


m 2 — mi 


uj { 0 A 1 N) (H { 0 A \(j)) ■,(/)). 


Clearly, 


,(A,JV) 

o, 


7 TJ 


(A) 


777-2 


7712 — 777i 
777-1 


7712 — 7711 


O 2 ); fa) - ^ 0 A + N) (#o A) (fa ); fa ) 

■ [ w 0 A to 7 t' ) (^0 A) ( < / , 2); ^2) - uj^ a ’*\h (a) ( fa)-, fa)] 

fa) - u { 0 A 1 N \h (A) (fa); fa)]. 


From the above argument, this right-hand side is exponentially small in the 
linear size L ^ of the system. Thus, in the sector of isolated eigenenergies 
from the rest of the spectrum with finite energy gaps, the twisted phase de¬ 
pendence of the arithmetic mean of the energies in the sector is exponentially 
small in the linear size 


B.4. Degeneracy of the sector of the ground state. By assuming that 
the degeneracy of the eigenenergies in the sector of the ground state is 
lifted when changing the value of the twisted phase </>, we will deduce a 
contradiction. In consequence, we prove that the degeneracy of the sector of 
the ground state is not lifted for any value of the twisted phase cj) E [0, 2tt). In 
other words, the sector of the ground state is (quasi)degenerate irrespective 
of the twisted phase (j). Clearly, by combining this with the result (IB.411) in 
Sec. IB.21 the proof of Theorem [2] is completed. 

Consider again the same situation as in Sec. IB.21 To begin with, we note 
that all of the eigenenergies are a continuous function of the twisted phase cj). 
Relying on this fact, we can assume that there exist a sequence of lattices, 
Ai C A 2 C ■ ■ •, and a sufficiently small phase cj)' E (0, 27r) such that for all 
the lattices A n , the phase cj)' satisfies the following two conditions: 

5E(fa) = max {5E(fa} 

*e( o,<yp 

and 

0 < 5E min < 5E(fa) < ^AE(0). 

Here, the phase 0' may depend on the sizes |A n | of the lattices A n , and 
the positive constant 5E mm is independent of the sizes |A n | of the lattices. 
Clearly, the second condition implies that the degeneracy of the sector of 
the ground state is lifted at the twisted phase <j) = <jJ. 

First, fix 0 = fa. Then, there exists a subsequence of lattices {A^} n C 
{A n } n such that the sector of the (/-fold ground state splits into i sectors each 
of which is c/j-fold degenerate for i = 1,2, and that there exist a nonva¬ 
nishing spectral gap 7 m between two adjacent sectors for m = 1 , 2 , — 1 . 

Clearly, the degeneracy and the gaps must satisfy 

l e-i 1 

'}2 ( h = q and ^2 7m < -A£(0), 

2=1 m= 1 
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respectively. We denote by Pq^’ N \(/)) the projection onto the z-th sector, 
and define 


u. 


(A,TV) 

0,i 


(••• ; ft) := ~Tr (• • • * = 1;2 


Qi 


Then, the above conditions about the spectral gaps are written 
(B.45) 


u ot5 ( H o A) (ft')* ft) - w$; N) (h^ a) (ft)-, ft) > n > o, * = i, 2,... 


I- 1 . 


Consider changing the twisted phase (j) from ft to ft\ € [0, <^> / ). Then, the 
deviation of the energy from the averaged energy for the z-th sector is given 

by 


AE^Oi) := ^ax {|S^._ i+j (^) -4 A ’^ ) (-ff l S A) (^l);0i)|}> 


* = 1 , 2 , 


where mo = 0 and 

i—1 

nii-i = qj for i = 2 , 3, ..., l. 

3 = 1 

We can find ft x € [0, ft) so that the spectral gaps between two adjacent 
sectors remain open when the twisted phase </> continuously varies from ft 
to (ft € [0, ft). Further, for a sufficiently small deviation ft — ft x of the 
twisted phase, we can assume that the deviation of the energy in the z-th 
sector satisfies 

(B.46) ^Eoftftft) < J 7j 

where 7 := min{ 7 i, 72 ,..., 7 £_ 1 }- Unless these spectral gaps close for vary¬ 
ing the twisted phase, the deviation of the averaged energy is exponentially 
small in the linear size L U) of the lattice as 

4 A ’ JV) (- ff o A) ( < ? !,, i); < / , l) -4^ ) (^o A) (^ , ) ;< ^) = 0( ex P[-Const.L (1) ]) 

which is the result in the preceding Sec. IB.31 Therefore, by combining 
this with the bounds (IB.451) and (IB.46p . one notices the following: If the 
condition (IB.461) is always retained during continuously varying the twisted 
phase from (j) = ft to (j) = 0 for any fixed lattice An, then the degeneracy 
of the low-energy sector must be lifted at ^> = 0. This contradicts the 
assumption that the energies of the low-energy sector are (quasi) degenerate 
at cj) = 0. Thus, there must exists a phase ft 2 £ [0, ft) such that the condition 
(IB.461) does not hold at the phase 4> = ft 2 . 

From these observations, we can assume that there exist a subsequence 
{An^}n °f th e sequence {A ^' 1 } n of the lattices, and the minimum value of 
the phase (j)\ € (0, ft) such that the condition (IB. 461) holds for all the sectors 
and that, in particular, there exists a sector A 1 ) that the deviation of the 
energy in the z^b-th sector satisfies 

AE o,im(fti) = 


(B.47) 
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for all the lattices A„ . Actually, since the eigenenergies are a continuous 
function of the twisted phase <p for a fixed lattice, there exists <p\ £ (0, <pf) 
for each lattice A® such that 

max{ A A (M (^i)} = 77- 

l<i<Z 4 

Therefore, there exists a subsequence {A^} n of the sequence {A,^°' ) } n that 
the condition (IB.471) holds for a sector . 

The condition (IB. 471) implies that the degeneracy of the i^-th sector 
is lifted at the twisted phase (p = <p\. Appropriately choosing the subse¬ 
quence {An ^} n of the sequence {A,^} n of the lattices, the i^-th sector 
splits into £{i^) sectors each of which is -fold (quasi)degenerate for 
i = 1,2 ,... and the spectral gaps 7^1) m between two adjacent sec¬ 

tors appear for m = 1,2,... , £(i^) — 1. Therefore, we can define the spec¬ 
tral projection P^’(^\{(p) and the expectation value w q^(i)^(' ■■;</>) for the 

(i^\i)-th sector for z = 1,2,... in the same way as in the above. If 

another sector j'W 7^ iW satisfies the same condition AE 0 ■ (i)O'i) = 7/ 4 for 

an infinite number of the lattices in the sequence {A^} n , then we make the 
same procedure as the i^-th sector. 

Further decrease the phase to <p = <p' 2 satisfying 0 < (p 2 < 4>\ < <p'. Then, 
we can define the deviation AE 0 70 ^(<^>2) °f the energy for the (*W,i)-th 
sector in the same way, and AE 0 -(p j((p 2 ) for the th sector. For a 

sufficiently small deviation cp^ — (p 2 , we can assume that 

(B.48) AE 0 , k (<p' 2 ) < ^7 

for the k -th sector satisfying 

AEo,k((p'i) < ^7 

and 

(B.49) A E q . w .{<P' 2 ) < | 7 (J (1) ) 

for the (j 4 -*, j)-th sector satisfying 

A^-cu^'i) = |7> 

where 

7(4 (1) ) == m in{7,-(i) ll> 7j(i) )2 , • • ■ ’ 7j(i)^(j(i))_i,7/ 4 } 

If these conditions are always retained during continuously varying the 
twisted phase from cp = cp^ to (p = 0 for any fixed lattices A„\ then the 
degeneracy of the low-energy sector must be degenerate at (p = 0. This 
contradicts the assumption again. 

If there exist a subsequence of the lattices and (p' 2 such that the condition 

AE 0tk (<p 2 ) = ^7 
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holds for a sector k in the first case ()B.48I) , then we make the same procedure 
as in the above zW-th sector. 

In the second case ()B.49[) . if there exist a subsequence {A^ } n of the 
sequence {An'*} n of the lattices and the twisted phase (j)' 2 such that the 
condition 

AE o = ^7(j (1) ) 

holds for the , j)-th sector, this condition implies that the degeneracy 
of the sector is lifted again. Therefore, we can make the same procedure. 

Since the initial degeneracy q is finite, the process must stop in finite 
times. Besides, once a spectral gap appears in the spectrum, it never closes 
because the total deviation of energy is bounded by 


1 , 1 , 

4 7 + 4^ 7 + 


1 , 1 , 

= 3 7 <2 7 ’ 


where r y' is the initial gap when the degeneracy of the sector is lifted. Thus, 
the degeneracy of the spectrum at the twisted phase = 0 is not recovered 
when continuously varying the twisted phase from (j> = <fi' to <J> = 0. This is 
a contradiction. Theorem [2] has been proved. 


Appendix C. Proof of the inequality (IB.101) 

Consider a Hamiltonian LfA) on the lattice A. We denote by the 
restriction of the Hamiltonian H^ onto the subset fl of the lattice A. Let 
A be a local observable, and define the time evolutions as 

A W(t) := e itHW Ae~ itHiA) 

and 

A W(t) := e itHm Ae- itHW . 

Then, one has [22] 

(C.l) \\A^ A \t) - A^\t)\\ < sgnt [ ds ||[(F (A) -i/ (Q) ),A( n) (f-s)]||. 

Jo 

In order to prove this inequality, we use the following identity: 

A^(t) - A^ n \t) = f ds ^-e isIiW A^\t - s)e~ isH(A) . 

Jo d s 

The derivative of the integrand in the right-hand side is computed as 

d JsHW 4(H)/. n -isH^ _ d j sH W J(t-s)H^ n ') A -i(t-s)H^ -isH^ 

-e A (t — s)e --e Ae 

= ie isH(A) [(H^ - H^),A^\t - s)]e- isH(A \ 


Combining these, the desired inequality (1C.II) is obtained. 
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Appendix D. Proof of the inequality (IB .261) 

The function f(z) in the integrand in the left-hand side of (1B.26|) is written 

f(T 2 e w ) = - 

77i=l n> 1 
xexp 

where we have written 


A(<j>) = a*P^ N \^)a 

for short. Using iT 2 e ie = iT 2 cos 9 — T 2 sin 0, the Schwarz inequality and the 
assumption on the spectral gap, one has 


f(T 2 e ie ) < ||af ||£(</OI|exp[-T 2 A£ mm sin0]. 


Therefore, the integral is evaluated as 
r f{z) 


-dz 


z — is 


< 


dd\f(T 2 e ie ) 

i wmw r 

■Jo 


T 2 


< a 


y/T 2 cos 2 9 + (T 2 sin 9 — s) 2 
d9 exp[-T 2 A£ min sin 9} 


V 1 - 2s / T 2 


<l|a|| 2 |TO)ll 


7r 


1 - e- T * AEU 


A E* 


T 2 


v'l - 2s/ T 2 ' 


This is the desired bound ()B.26D . 


Appendix E. Proof of the inequality (IB.301) 

Note that 
(E.l) 

772=1 ? 2>1 

x exp [i(E§2.(<l>) - E^ N \(j)))t], 

where we have written 

A\4>) = aP^’ N \cj))a* 

for short. Therefore, the integral with respect to time t in (IB. 271) is written 

(E.2) — it — 

Zm J_ T2 t — is 

In order to estimate this integral, we prepare a tool 12 0123]. 

For Eel, one has 

i r+° 

JEt -at 2 _ 1 / 

2 V^J-oo 


dw e lwt exp[— (w — E) 2 /{ 4a)]. 
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Therefore, 


(E.3) 


r+Tz giEtg—at 2 


2ni 


—T 2 


t — z 


dt = 


1 


f +00 


2^[Ea J_ 

1 

27ri 


dw exp[— (w — E) 2 /(4a)] 


r+T 2 „iwt 

/ dt- - 

I—T 2 t ^ 


for z E C satisfying Im z > 0. For the integral with respect to time t, we 
write 


(E.4) 


= hj[ 


+T 2 p iwt 

dt- - 

-T 2 t-z 


with Im z > 0. Then, the integral (IE.3D is written 
(E.5) 


_L [ 

2iri J_ T2 


+T 2 at 2 -y 


t — Z 


-dt = 


2y/ira J _ 

1 


/ +00 

dw exp[— (w — E) 2 /(4a)\F Zj T 2 (w) 

-OO 


L 


OO 

+00 


2 y/Wa. 

+ R\ + R‘2 


dwe iwz ex.p[-(w-E) 2 /(4a)\ 


with 


and 


Ri = 


1 r +00 

0 7 = dw exp[-(u> - £) 2 /( 4 «)] [F Z) t 2 {w) ~ e lwz ] 

2^/txol J 0 


R2 = 


1 


2yJiroi 7„ c 
One can easily obtain 


dw ex.p[—(w — E) 2 /(4a)]F z> T 2 (w). 


(E. 6 ) 

and 


1 


\F z T2 {w) - e iwz \ < [1 - e~ wT2 ] , for w > 0 

wT 2 


\ f z,t 2 { w )\ < 


1 


l _ g-MT 2 


, for w < 0 , 


|w|T 2 

where we have assumed Im z > 0. Using these bounds, one has 

r»+oo 


1 /*+°° 1 
l^il < 77 = / dw exp[-(w-E) 2 /(4a)]—- [l 
2y/Tta J 0 


— e 


—wT 2 1 


and 


l^|< 


2^/twl J_ 


f 


dw exp[— (w — E) 2 /(4a)]- 


w\T 2 


l _ e ~\ w \ T 2 


Further, applying the Schwarz inequality yields 

r+oo 


1 /*+oo 1 r 

\Ri\ + I-R 2 I < 77 = / dw exp[-(w - E) 2 / (4a)]- — 1 - e “ HT2 
AyJ'KOi J_ OG 17/71 _Z 2 


< Const. 


aV 4 T 2 1/2 ' 


Thus, Ri and R 2 are vanishing in the limit T 2 —>• 00 . 
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From (IE. ID and ()E.5D . the integral I\ is written 


h = 


2iri 


r+T 2 


'-t 2 


dt - —uji A,N \B^ A \t;(j))A\(f));(l))e at 


t — is 




m= 1 n> 1 
^ P+T2 


2ni 


/ \1 2 1 

T dt - 42(^))] e_Qt2 

= ^ E E (^ff(^)^(^)^ ) (0))(^ ) (0)^'(^)^S(0)} 

^ ra=l n>l 
i /*oo 

x ^7= / ^e“" ;s exp[-{u; + (E;( l iV) ((/.)-Sj^( ( /)))} 2 /(4a)] 
lAV 71 ^ Jo 

+ i?i + i ?2 

Therefore, we have 

161 < i||BWIIM'W)l| £ - <aE " i " )1/,to) < ^|BWIII|a|| 2 e-‘ aE "'"» a /(-«) 

in the limit T 2 —>• 00 . Here, we have used the assumption ()B.21D on the 
spectral gap. 

Appendix F. Proof of the inequality (IB.31D 

Let 5 be a large positive number which satisfies S < T 2 . We decompose 
the integral I 2 into three parts as 


h = 


1 

27 xi 


r+T 2 


—T 2 


dt-±-4 A ’ N \a*p( A ’ N \ct>)[a, B^(t- J>)\A)e- at2 


t — is 


— !•) + ^2 + I 2 


with 


/o' = 


2ni 


—s 


-t 2 


dt—^—ul A ' N) (a*P 0 (A,7V) 0) [a, P (A) ( t ; <f >)]; 4>)e 


-at z 


t — is 


I°2 = / + %t^4 A ’ JV) (a^o (A ’ JV) (^)[«^( A )(t^)];^)e-“ t2 


27TZ 


and 


r-+S 

di 

i-s t — is 

r+T 2 


It = [ 2 dt-^-u^^P^ma, B^(f, </>)]; 

ZtlTlj J15 £ ZS 


The hrst and third integral are estimated as 


\I?\ < 


1 1 


|af||P(0)||e-“ s < 




D - a s 2 


2 27t ^52 _|_ s 2 y a 11 11 11 vv ' ;n 2y/Wa S 

For the second integral /l,, we use the Lieb-Robinson bounds |3[23|. Note 
that the operator />(</>) is a sum of the local operators bj((f>) as 

£(</>) = Ew 
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From this and the Lieb-Robinson bounds, one has 

||[a,R (A) (i;</>)]|| < Yl\\ 0)]|| 

3 

< Const.|supp a|||a|| E |supp bj(4>)\\\bj{4>)\\e ^ r (e vW - 1) 

3 

< Const.|supp a|||a||||R(</>)||e _Al? ’(e 1 '^l — 1), 

where the operator (i; J>) is the time evolution of bj (0). Using this bound, 
the integral I® can be estimated as 

I/ 2 I < Const.|supp a|||a|| 2 ||-B(</>)||e 


2\\nfA\W--nr I df 
>-S 


+S e v \ t \ _ Y 


,-ar 


r+S vt _ 1 

< Const.|supp a| ||a || 2 ||R(^)||e _Mr / dt - 

Jo t 

< Const.|supp a| ||a|| 2 ||.B(^>)|| exp[u5 — fir]. 


Combining this with the bounds for 1^. one obtains 


|/2 |<||a|| 2 ||R(</>)|| 


1 e 


-«S 2 


/7r a S 


+ Const.|supp a| exp [u S' — fir] 


Appendix G. Proof of the inequality (IB. 32(1 
Note that 

(G.l) f(t) = uj ( 0 A ' N) (A^B^ (t; (j))- J)) 

=- it 


m = 1 n > 1 


x exp [i{E^ N) ( 4 >) - 

Therefore, the integral with respect to time t in I 3 is 


(G.2) 


1 r +T 2 1 22 

271 j _ T2 d t— s eM^E n , m t\(e<" - ), 


where we have written AE nm = E^\<J>) — (r/>) for short. From (IE. 41) 

and (IE. 61) . one has 

1 f +T2 1 2 

— dt - -— exp[iAE nrn t]e as 

2m J-t 2 t — is 

= exp [as 2 - AE n:7n s\ + R :i , 

where the correction R 3 is vanishing in the limit T 2 —>• 00 and we have used 
the assumption (IB.2111 on the spectral gap. Since 


1 


2^/na 


/ +00 

e~ sw exp[—(u; - AE n>m ) 2 /(4a)] = exp[as 2 - AE n , m s], 

-OO 
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the above integral is written 


1 f +T2 , 1 

f at -— ex 

—T2 t — hs 


2 tti 


\)[iAE nrn t]e 


Y r+o o 


2 ^pna J_ 

Combining this with (IE.51) . one obtains 


e s1 "exp[-0 - AE n . :m ) 2 /(4a)\ + R 3 . 


1 

2 iri 


r + T 2 1 „ 

/ dt - — r "- Ar + u ~ as ~~ at 

l_T 2 t — is 

^ r +00 


■ exp[iAE n>m t](e“ s2 - e “* 2 ) 


2y/ira J_ 

1 


OO 

■+OO 


2 y/Tra 


/ +00 

e~ sw exp[—(w; - AJ?„ im ) 2 /(4a)] + R 3 

-OO 

dw e~ sw exp [—(in - A£' nim ) 2 /(4a)] - R 3 - R 2 

1 r° 

= 0 7 — / e _, 5 U ’exp[—(in — A£ , ri)rn ) 2 /(4a)] — R\ — R 2 + R 3 - 

Ay/ircx J —00 

Using 0 < s < cr and AE n , m > AE mm , and choosing a = AF min /(2c?’), the 
above integral in the last line can be evaluated as 

—U [° e~ sw exp[—( W -AE n , m ) 2 /(4a)]<i e - AsminCT / 2 . 

2 \/vra J-00 2 

Therefore, from (1G.1I) and (IG.2D . the integral I 3 of (IB. 291) is estimated as 


i/ 3 i < \hf\\Bme - AEm '^ 12 


in the limit T 2 —>• 00 . 


Appendix H. Differentiability of 0( 2 )(</>i) 

Fix cj )® E [0, 2ir), and consider cf) 1 in the neighborhood of </>j°\ Then, 
the product C^ 2 \(j)i)C^ of the two unitary matrices is nearly equal 
to the identity matrix. Consider the Hermitian matrix A0( 2 '(0i) which is 
defined by 

exp [iA0( 2 )(0i)] = C^ 2 \(f)\)C^ 2 \ 4 ^ > )^ ■ 

For a sufficiently small \(f>\ — 7>^|, the matrix A0( 2 )(7>i) is well defined and 
differentiable with respect to (f>\. Clearly, one has 

(H.l) 1 ) = exp [iA0®(^)] C (2) (4 0) ). 

Further, one obtains 

(H - 2 ) i 

C ’ ^ C^A-C^) = < jfJ ’ Tr Aa®«(*) 

= i[Tr A0 (2) (^ 1} ) - TrA0 (2) (^ O) ) ] 

= *Tr A0®(^ 1} ). 
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On the other hand, using the expression 

C (2) (<M = exp [i0 (2) Oi)] 

and taking the determinant for both sides of (IH.ll) . we have 

Tr 0 (2) (</>i) = Tr A0^(0i) + TV 0 (2 )(^ O) ) (mod 2 tt). 

This implies 

TV A0^(^ 1} ) = Tr 0( 2 )(<^ 1} ) — Tr 0®(0^ O) ) (mod 2vr). 
Substituting this into the right-hand side of (IH.2D . we obtain 

/t #iTrC( 2 )(^)tAc( 2 )(^) 

J 4 0) 

= ?:[TV 0 (2) (^ 1} ) -Tr 0 (2) (^ O) )] (mod 2vr). 
This result justifies the formula (16.91) . 


Appendix I. Proof of Lemma H] 

To begin with, we recall that the relation between the two projections is 
given by 

P^ A ' N) (if), a) = e^[-ia<l)in x ]P^ A ’ N \4>)exp[ia(l)in x ]. 

By differentiating with respect to 0i, one obtains 

9 = exp[-ia(j)in x ]P^ A ’ N \&)exp[ia(j)in x ] 


901 


— ( ian x ) exp[—ia(j)in x \PQ A ’ N \(j)) exp[fa0in 


+ exp [~ia(j)in x \P^ A ' N \(j)) exp[iacj)in x \(ian x ). 
Using these relations, we have 

Tr P 0 (A,W) (0, a) [Po ( t’ W) ( 0»,«), -f 0 , 2 ’^(0> “)] 

= Tr p( A ' N \fr[P$ N \4>),P$ N \j>)} 

- ia Tr P^ N \*)[n^ N \<t>),P^\4>)] 

+ ia Tr P o (A ’ Ar) (0)[P o (A ’ JV) (0)u, x , P^ ,N \^)}. 
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Therefore, it is enough to handle the second and the third terms in the 
right-hand side. Except for the factor ia, these two terms are written 

Tr p( A ’ N \j>)P$ N) (j>)n x - TV p( A ’ N \&n x P^ N) (j>)P$ N \j>) 

+ Tr A ’ N \i) - Tr P^’ N) (^)P^ N) (^)P^ N) (^)n x 

= Tr p( A ’ N \j>)P$ N \j>)[ 1 - 
+ Tr P^ N \^)n x [ 1 - P( A ’ N \&]P$ N \& 

= E [<^S(^)^o ( , 2 ,7V) (^)[1 - i- o (A ’ 7V) (0)]^ci>^(^)> 

m= 1 

+ (dw.M 1 - d A, '’ ) («]<2 , ' v) (i)-i>Sw>]. 

where we have used Proposition [31 Further, by using (16.41) . this is written 


E [! - P 0 A ' N) (<t>)\ n x® { 0 N m(<t>) 

m= 1 


E E *©(*))<*©(*■ ^*S2(« 

m=lm' = l 2 

= E e-{Cw.".*£®) 

m=l ^ 

- E E 


m 

+ 

Q 


Clearly, the double sum in the right-hand side is vanishing. Integrating the 
first sum with respect to 4>2 yields 


/» 27t 9 a 

J o E 

= E [(^£2(^i. 27r )> n x^oS(^i> 27r )) - ($K(0i>°)> n *$S2(0i«°)) 

ra=l 
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By using the relation (16.81) . we can show that the right-hand side is vanishing. 
Actually, we obtain 


( 0 ^> 2 ^ 0 ^ 27r )) 

m= 1 

= £ £ £ 

m=l m'=l m"=l 



m'=l 


Thus, the corresponding contribution is vanishing. 
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